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ABSTRACT 



The advent of modern infrared astronomy has brought into evidence the role 
played by the interstellar dust in galaxy formation and evolution. Therefore, to fully 
exploit modern data, realistic spectrophotometric models of galaxies must include this 
important component of the interstellar medium (ISM). 

In this paper, the first of a series of two devoted to modelling the spectra of galaxies 
of different morphological type in presence of dust, we present our description of the 
dust both in the diffuse ISM and the molecular clouds (MCs). 

Our galaxy model contains three interacting components: the diffuse ISM, made 
of gas and dust, the large complexes of MCs in which active star formation occurs and, 
finally, the populations of stars that are no longer embedded in the dusty environment 
of their parental MCs. 

Our model for the dust takes into account three components, i.e. graphite, silicates 
and polycyclic aromatic hydrocarbons (PAHs). We consider and adapt to our aims two 
prescriptions for the size distribution of the dust grains and two models for the emission 
of the dusty ISM. We cross-check the emission and extinction models of the ISM by 
calculating the extinction curves and the emission for the typical environments of the 
Milky Way (MW) and the Large and Small Magellanic Clouds (LMC and SMC) and 
by comparing the results with the observational data. The fin al model we have adopted 
is an h ybrid one which stems from combining the an alysis oflGuhathakurta &: Draind 
l(T989|) for the emission of graphite and silic ates a nd Pugc t et alJ lll985|) for the PAH 
emission, and using the dist ribution law of Wcingar tnerfc Draind l)2001a|) and the 
ionization model for PAHs of lWeingartner fcDrainel l|2001b|h 

We apply the model to calculate the spectral energy distribution (SED) of Single 
Stellar Populations (SSPs) of different age and chemical composition, which may be 
severely affected by dust at least in two types of star: the young, massive stars while 
they are still embedded in their parental MCs and the intermedi ate- and low- mass 
AGB stars when they form their own dust shell around fsee lPiovan et al1l2003l for 
more details about AGB stars). 

We use the "Ray Tracing" method to solve the problem of radiative transfer and 
to calculate extended libraries of SSP SEDs. Particular care is paid to model the 
contribution from PAHs, introducing different abundances of C in the population of 
very small carbonaceous grains (VSGs) and different ionization states in PAHs. The 
SEDs of young SSPs are then compared with observational data of star forming regions 
of four local galaxies successfully reproducing their SEDs from the UV-optical regions 
to the mid and far infrared (MIR and FIR, respectively). 

Key words: ISM: dust, extinction - infrared: ISM - galaxies: Magellanic Clouds, 
radiative transfer. 



1 INTRODUCTION 

The interstellar dust, either tightly associated to stars 
and/or dispersed in the ISM, has got more and more atten- 
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tion over the years, in particular with the advent of mod- 
ern infrared satellites (e.g. IRAS, COBE and ISO), because 
of its role in m any astrophysical phenomena (see iDraind 
l2003allbHl2nol for more details). 

Leaving exceptions aside, there are at least three main 
circumstances in which dust influences the stellar light: (i) 
It is long known that for a certain fraction of their life very 
young stars are embedded in the parental MCs. Even if the 
duration of this obscured period is short, its effect on the 
light emitted by these stars cannot be neglected as a signifi- 
cant fraction of the light (initially almost all) is shifted to the 
IR region of the spectrum, (ii) Low and intermediate mass 
stars in the Asymptotic Giant Branch (AGB) phase may 
form an outer dust-rich shell of material obscuring and re- 
processing the radiation emitted by the star underneath, (iii) 
Finally, thanks to the contribution of metal-rich material by 
supernovae and stellar winds in Wolf-Rayet and AGB stars, 
the ISM acquires over the years a dust rich component. The 
UV-optical radiation emitted by stars passing through this 
dust-rich intergalactic gas is absorbed and then re-emitted 
in the far IR. 

In this paper, first we develop a model for the ab- 
sortiom/emission properties of a dusty medium and second 
we apply it to derive the SED of young stars s till embedded 
in th eir parental MCs. In a companion paper jPiovan et alJ 
2005), we will present detailed chemo-spectrophotometric 
models of galaxies of different morphological type whose 
SED from the UV to the far IR is derived including dusty 
MCs and the presence of diffuse, dust-rich ISM. 

Stars are preferentially born inside massive, dense and 
cool MCs characterized by low temperatures (T < 9 — 15K), 
masses in the range from ~ 10 4 to 10 6 Mq and dimensions 
from ~ 6 to 6 0pc, the most massiv e and big ones being 
less numerous (ISolomon et al.lll987l) . Furthermore all re- 
gions with active star formation, e.g. in the MW, are also as- 
sociated to H2 clouds, thanks to which the MCs are mapped 
by means of radio surveys. In general they seem to be or- 
ganized in hierarchical structures forming very complicate 
and large complexes, inside which a large number of sub- 
structures of higher density are found. The observation of 
the earliest evolutionary phases of young stars is severely 
hampered by the presence of dust, which absorbs and dif- 
fuses a large fraction of the radiation emitted by the stars in 
the optical returning it in the MIR/FIR. However, as soon 
as the first-born, massive stars evolve, their strong stellar 
winds and mass-loss, intense ionizing radiation fields and fi- 
nal explosion as type II supernovae will eventually destroy 
the MCs in which they are embedded. Massive stars become 
eventually visible in UV-optical regions of the spectrum. The 
time scale for this to occur is ~ 10 6 — 10 8 years, i.e. the life- 
time of the most massive stars in the population. To summa- 
rize, we are able to map the location of MCs thanks to radio 
data, to follow the first evolutionary stages of star forma- 
tion thanks to FIR data and finally to directly observe the 
stars in the UV-optical when the parental clouds have been 
evaporated. However all intermediate stages are precluded 
to direct observations because they are obscured by dust. 
They must be inferred from theory, which unfortunately is 
still far from being fully satisfactory. 

The problem is particularly severe when the evolu- 
tionary population synthesis technique (EPS) is applied 
to model the SED, the integrated magnitudes and broad- 



band colors of galaxies, by folding the properties of SSPs 
of different age and chemical composition on the star for- 
mation history (SFH). Classical studies of this subject 
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SSPs. In other words, "bare" spectra of SSPs are used to 
synthesize the SED of a galaxy. This approximation sounds 
acceptable when modelling the SED of old systems, such as 
early-type galaxies, in which star formation took place long 
ago, even if also in this case t he presence of dust around AGB 
stars should not be ignored JPiovan et alJl2003t: iTemi et alJ 
2004). This is certainly not the case with the late- type and 
starburst galaxies that are rich of gas, dust and stars of any 
age. Therefore, it is mandatory to include the effect of dust 
on the radiation emitted by young stars. Because of the high 
density in the regions of star formation, the optical depth 
may be very high also for IR photons and the full problem 
of radiative transfer has to be considered. 

The plan of the paper is as follows. In Sect. [5] we 
model in detail the extinction and emission properties of the 
dusty ISM. In particular we present the optical properties of 
graphite, silicates and PAHs, the two distribution laws for 
the grain sizes (shortly named MRN ad WEI models), the 
cross sections and the dust to gas ratio. In Sect.|3]we address 
the topic of emission from dust grains, and present two mod- 
els shortly indicated as GDP and LID. The theory is applied 
to derive the average extinction curves and emission proper- 
ties for the diffuse ISM of the MW, LMC and SMC ( Sects. H 
14.11 and 14.21 respectively). Then we apply our best solution 
for the ISM extinction and emission to derive the SEDs of 
young dusty SSPs. We summarize the fundamentals of spec- 
tral synthesis in Sect. Eland in Sect.HJwe present our choice 
for the spatial distribution of young stars inside the MCs. In 
Sect.|7|we cast the problem of radiative transfer introducing 
the key parameters and discussing the dependence of the op- 
tical depth on the cloud physical parameters. In addition to 
this, we quickly summarize the the "ray tracing" technique 
to solve the problem of radiative transfer throughout an ab- 
sorbing/emitting medium. In Sect.|H]we present the SEDs of 
very young SSPs surrounded by their MCs. Having in mind 
the application of our results to studi es of population syn- 
thesis in galaxies dPiovan et al-l feoOS). we have calculated 
a library of young SEDs with dust for large ranges of the 
parameters, paying particular attention to the emission of 
PAHs. In addition to this, we show results for a simple way 
of modelling the evaporation of the MCs surrounding young 
SSPs. In Sect. |U] we compare our SEDs with observational 
data of star forming regions. To this aim, we considered the 
central regions of four starburst galaxies, namely Arp220, 
NGC253, M82 and NGC1808. Finally, in Sect.EHlwe draw 
some general remarks and conclusions. 



2 MODELLING THE PROPERTIES OF DUST 

Chief workbenches for the study of the physical properties 
of the interstellar grains are the extinction curves and the 
IR emission of dust observed with growing detail in different 
physical environments. From the characteristic broad bump 
of the extinction curves in the UV a t 2175 A and the ab- 
sorption features at 9.7/im and 18^m forainel2003J) we can 
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already infer that a two components model made of graphite 
and silicates is required. The study of the emission adds fur- 
ther constraints. A population of VSGs has been invoked to 
reproduce the emissions observed by IRAS in the pass-bands 
12/j,m and 25p,m. The VSGs temperatures can fluctuate well 
above 20 K if their energy content is small compared to 
the energy of the a bsorbed phot o ns llLeeer fc Pueetlll984t 
Desert et al.1 I198& Ipwekl Il986l: iGuhathakurta fc Drainel 



19891: I S^b^nrn^rTOri^ri^ruegeir 119921 : IDraine fc Lil 1200 



Li fc Drainel 1200 if) . Excluding that VSGs are made of sil- 



icates simply because the 10/xm emis sion feature of sili- 
cates is not detected in diffuse clouds jMattila et al.lll996t 
lOnaka et aHll996Tl . most likely they are made of carbona- 
ceous material with broad r anges of shapes, dimensions, and 
chemic al structures (see also lDesert et alJl986HLi fc Drainel 
2002a, for a discussion of this topic). 

Emission lines at 3.3, 6.2, 7.7, 8.6, and 11.3/wn, origi- 
nally named unidentified infrared bands (UIBs), have been 
firstly observed in lu minous reflection nebulae, planetar y HII 
regions and nebulae dSellgren et al 1 ll983tlMathis!ll99rJ) and 
subsequently also in the di ffuse ISM wi th IRTS ([Onak a et alJ 
119961: banaka et ail ll996Tl and ISO llMattila et alJll996Tl . 
There is nowadays the general consensus that these lines 
owe their origin to the presence of PAH molecules, vibra- 
tionally excited by the absorption of a U V-optical photon 
jLeeer fc Pueetl Il984t iLi fc Drainel 1200 J) . Currently these 
spectral features are commonly referred to as the aromatic 
IR bands (AIBs). 

Based on these considerations, any realistic model of a 
dusty ISM, to be able to explain the UV-optical extinction 
and the IR emission of galaxies, has to include at least three 
components, i.e. graphite, silicates, and PAHs. Furthermore, 
while it may treat the big grains as in thermal equilibrium 
with the radiation field, it has to allow the VSGs to have 
temperatures above the mean equilibrium value. In order 
to obtain the properties of a mixture of grains, we have to 
specify their cross sections, their dimensions, and the kind 
of interaction with the local radiation field. 

Optical properties. The optical properties of PAHs, 
silicate and graphite grains together with the corre- 
sponding dimensionless scattering and absorption coeffi- 
cients, Qsca(a , X) a n d Q a i, B (a, A) ha v e bee n taken from 
IDraine fc Led Jl984Tl. iLaor fc Drainel dl993T) . IDraine fc Lil 
fcOOlT l and lLi fc Drainel fcOOll) . These coefficients are de- 
fined as the ratio of the cross section a to the geometrical 
section ira 2 , where a is the dimension of the grain. 

Distribution laws: MRN-like and WEI models. 

In the ISM the grain dimensions are likely to vary over a 
large range of values. Therefore in order to properly model 
the optical properties of the ISM we need to specify the 
distribution law of the grain dimensions and to fix their 
upper and lower limits. Two cases are considered. 

In the first one, shortly refer red to as "WEI", we 
adopt t he analytical law proposed bv lWeingartner fc Drainel 
feOOlaT) . It is a very complicate relationship for the size 
distribution of carbonaceous dust grains, which simultane- 
ously deals with graphite and PAHs and smoothly shifts 
from PAHs to small graphite grains, considering PAHs as 
the extension of carbonaceous grains down to the molecular 
regime: 



1 dn a „ . . C / a \ ag „ , 

— = D a +— F(a,/3 g ,a t , g )x 

n H da a \a t , g ) 
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exp {— [(a — a t , g ) /a c , g \ ) a > at, g 

where C g is the abundance of carbon, a±, g is a transition 
radius securing a smooth cut-off for dimensions a > at, g , 
parameter that controls the cut-off steepness, and 
a g is the exponent of the power law. D (a), the sum of two 
log-normal terms that helps to better reproduce the emission 
by very small carbonaceous grains, is given by 
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for a > 3.5 A. Following I Weineartner fc Drainel (^OPlaT) . the 
term B; is 
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3 exp (-4.5 a 



(2tt 
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where me is the mass of a C atom, p — 2.24g cm is the 
density of graphite, bc,i = 0.75 be, be, 2 = 0.25 be with be 
the total C abundance (per H nucleus) in the log-normal 
populations, ao,i = 3.5 A, 00,2 = 30 A and a — 0.4. The 
function F (a, P g ,at, g ) in eqn.^is a suitable correcting term 
of curvature [see eqn. (6) in IWeingartner fc Drainel <l2001aTl 
for more details]. Finally, to obtain the total abundance of 
carbon we need to add the contributions coming from the 
two log- normal distributions, i.e. D(a), to the term con- 
trolled by C g . An interesting property of eqn. Q is that 
thanks to its ad-hoc analytical form very good fits of the 
extinction curves are possible by varying the contribution 
of the very small carbonaceous grains to the C abundance 
in the ISM, a quantity not yet firmly assessed. As shown 
bv lWeingartner fc Drainel <l2001aT) . the sole extinction curve 
cannot constrain the carbon abundance be, but only provide 
an upper limit. For any ratio of visual extinction to redden- 
ing Rv = A (V) / E (B — V), the limit is reached when the 
very small carbonaceous grains (PAHs and small graphite 
grains) are able to account for all the ultraviolet bump of 
the extincti on curve. A relationsh i p simil ar to eqn. Q is 
adopted in IWeingartner fc Drainel (l2001aT) for the silicate 
grains but, since there is no need for log-normal distribu- 
tions of VSGs, the abundance of silicates is completely de- 
scribed by the parameter C„. With the a id of the above 
relations bv lWeingartner fc Drainel i2001aT) it is possible to 
reproduce the average extinction curves of the diffuse and 
dense ISM of the MW, LMC and SMC. 

In the second case, shortly referred to as "MRN", we 
have modifie d and extended the power-law proposed by 
iMathis et al.l lll977f) . To better explore all the possible situ- 
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ations, we split the distribution law of the i-th component 
in several intervals labelled by k 



1 drii 
tih da 
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with the condition that a 



P+Pi 



p+fy+t Pj+i where 



j = 1, ...,fc — 2 if k > 2. The meaning of a Cfc is obvious, 
whereas the a bu are the "m odulation factors" introduced by 
iDraine fc Andersonl il985Tl . For k = 2 relation @ reduces to 
a single power-law distribution with o cl = a m in and a Cfc = 
a-max- Usually, k = 2 for PAHs, and 2 ^ k ^ 4 for silicates 
and carbonaceous grains. In most cases we adopt k — 2 
or k = 3. More complicated power laws are obtained with 
k = 4. 

Cross sections The cross sections of scattering, ab- 
sorption and extinction for the grain mixture in each com- 
ponent of the ISM are defined as 



na Qp,i(a,\) da (5) 



da 



where the index p stands for absorption (abs), scattering 
(sea), total extinction (ext) and the index i identifies the 
type of grains, amin,i and a max ,i are the lower and up- 
per limits of the size distribution for the i-type of grain, 
and finally uh is the number density of H atoms. The 
total dimension- less extinction coefficient is Qext (a, A) = 
Q abs (a, A) + Q sca (a, A). 

With the aid of the above cross sections it is possible to 
calculate the optical depth r p (A) along a given path 



SMC JPeHl 992l. A simple way of incorporating the effect of 
metallicity is to assume that 8 cc Z in such a way to match 
the approximate average results for MW, LMC and SMC: 
5 = <5q (Z/Zq). This relation simply implies that metal-rich 
galaxies are also d ust-rich. The above relation 8 cc Z agrees 
with the results bv lDwekl <l 19981) based on evolutionary mod- 
els for the compositions and abundances of gas and dust in 
the MW. However, it oversimplifies the real situation. 

The slope of the extinction curves greatly varies pass- 
ing from the MW to LMC and SMC. Moreover, the bump 
due to graphite de creases from MW to LMC and SMC 
JCalzetti et alJl994l) . These differences have been attributed 
to the metallicity decreasing from Z = Zq in the solar vicin- 
ity, to Z = \Z Q in the LMC, and Z = \Z Q in the SMC. 
However, the metallicity difference does not only imply a 
difference in the absolute abundance of heavy elements in 
the dust, but also a difference in the relative proportions, 
i.e. in the composition pattern. Probably the ratio graphite 
to silicates varies from galaxy to galaxy. Despite these un- 
certainties toevriendt et alj | l999h . the relation 8 cc Z is 
adop ted to evaluate t he amount of dust in galaxy models 
(e.g. ISilva et alJll998T) by simply scaling the dust content 
adopted for the ISM of the MW to the metallicity under 
consideration. In this study, using the data for the extinc- 
tion curves and emission of the MW, LMC and SMC, we 
seek to build a model of the dusty ISM describing the effect 
of different metallicities in a more realistic way. It is clear, 
however, that the problem remains unsettled for metallici- 
ties higher than the solar one, for which we have no infor- 
mation. The relation 8 cc Z together with the relative pro- 
portions holding good for the MW extinction curve, can be 
reasonably adopted also for other dusty ISMs characterized 
by metallicities close to the one of the MW. 



t p w = a P' i w / nndl = a P' i w x Nh ( 6 ) 

where L is the optical path and the meaning of the other 
symbols is the same as before. In this expression for the vari- 
ous T P (A) we have implicitly assumed that the cross sections 
remain constant along the optical path. 

Dust to gas ratio The coefficients be, C g , C s for the 
WEI model and A, for the MRN model fix the abundances 
of carbon and silicates with respect to the abundance of hy- 
drogen. In order to apply our models to a wide range of 
metallicities and ages, it is important to discuss how these 
coefficients would change as function of the chemical compo- 
sition of the ISM (the metallicity at least). The dust-to-gas 
ratio is defined as 8 = Md/Mn where Md and Mh are the 
total dust and hydrogen mass, respectively. If m% is the to- 
tal mass of the i-th type of grains then Md = J^. rrii, where 
the mass rrn is obtained from integrating over the grain size 
distribution (which depends also on the abundance coeffi- 
cients). Knowing the gas content of a galaxy (in principle 
function of the age), the amount of dust in the ISM depends 
on 8. For the MW and other galaxies of the Local Group, 
8 is estimated to vary from 1/100 to 1/500. In models for 
the MW, 8 = 0.01 is often adopted, while 8 = 0.00288 and 
8 = 0.00184 are typical values for the LMC and SMC. The 
mass-density ratios of interstellar dust are roughly in the 
proportion 1:1/3:1/5 going from the MW to the LMC and 



3 EMISSION MODELS 

In the following, two cases are considered to evaluate 
the emission of a dusty ISM. First, the easy to calculate 
model that stems from the ther mal-continuous approach by 
iGuhathakurta fc Draind lll989f) for the temperature fluctu- 
ations of silicate and graphite gra i ns and a modification of 
the prescriptions by IPueet et alJ 1119851) for the tempera- 
ture fluctuations of PAHs. The calculation of the ioniza- 
tion state of PAHs is also p roperly introduced following 
IWeingartner fc Draind J2001al) . We will refer to it as the 
"GDP" model. Second, a more elaborated model in which 
we implement the state-of-the art for the temperature fluc- 
tuations of the ve ry s mall grains as recen tly proposed by 
IDraine fc Li (1200 j) and lLi fc Draind (1200 ll) . This model will 
be shortly indicated as the "LID" model. 

3.1 The GDP model 

Emission from silicate and graphite grains. In order 
to calculate the temperature distribution acquired by small 
spherical grains of dust heated by photons and collisions 
with energetic particles, it is useful to define a vector state 
P (t) . The state of a grain is determined by its enthalpy E 
and the components Pu (t) of the state vector are the prob- 
ability for a grain to find itself in the fc-th bin of enthalpy at 
the time t. The enthalpy spans the range [i? m i n , -Emax] which 
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in turn splits into TV discrete bins each of which character- 
ized by Ek and AEk, i.e. the average value of the enthalpy 
and the width of the entha lpy bin. The enthalpy grid we ha ve 
adopted is the same as in lOuhathakurta fc Draine) 0989). 

Let us also define the transition matrix Af t i represent- 
ing the probability for unit time that a grain in the state i 
may undergo a transition to the state / 7^ i. Neglecting the 
transitions to enthalpies outside the interval [£ m i n ,£ mMt ], 
the condition of statistical equilibrium (steady state) is given 
by dP f* /dt = jGuhathakurta fc Drainell989tlDraine fc Lil 
2001). Taking into account all the transitions to and from 

JV 

any bin we have ^2 Af t iPf 3 = 0. Together with the nor- 

<=i 

malization J^ i=1 Pi 3 = 1 condition, it forms a system of 
linear equations whose solution yields the equilibrium tem- 
perature distribution. All technical details on the solution of 
this steady state system and the transition ma trix Af.j can 
be fou nd in lGuhathakurta fc Drainel (ll98Sl) and lDraine fc Lil 
fcOOll) . 

The enthalpy and temperature intervals are usually cho- 
sen in such a way that the distribution function Pk (T) 
smoothly changes from one bin to another and the en- 
ergy balance betwee n absorption and emission is conserved. 
Lis fc Leunel Jl99ll) use a grid with TV ~ 200. However 
Siebenmorgen fc Kruegel il992t) note that grids with high 



TV do not necessarily improve the accuracy of Pk (T) / AEk 
and accordingly adopt TV — 60, even if TV ~ 30 could be 
suffici ent in some cases. In the latest models of lDraine fc Lil 
<200ll) the number TV of bins is in the range 300 < TV < 1500, 
each bin (the first ones in particular) contains at least two vi- 
brational states, the bin widths smoothly vary and the last 
bin is chosen at high enough energy so that the expected 
population is very small (probability Pn < 10~ 14 ). Here we 
adopt TV ~ 150 — 200, which seems to secure computational 
speed and accuracy at the same time. 

A question soon arises: below which dimensions dfi u 
must the temperatu re fluctuations be taken into ac- 
count? According to [G uhatha kurta fc Drainel dl989f) and 
ISiebenmorgen fc KruegeTi^92T) . the typical dimension sep- 
arating small from large grains is about cifiu ~ 150 A in- 
dependently of the radia tion fie l d. How ever, this conclusion 
was questioned bv lLi fc Drainel i200 ll) who found that also 
for dimensions of the order of afiu — 200A the tempera- 
ture distribution is sufficiently broad thus introducing some 
uncertainty in the estimate of the emission. To take this re- 
mark into account we adopt here afi u ~ 250A. Finally, the 
emission j^ mal1 for small graphite and silicate grains is given 
by 



■ small 



a flv 



a 2 Qabs(a,X)B x (T(a)) x 



x ?m dT ±^ da 



dT 



Tin da 



(7) 



where dP (a) /dT is the distribution temperature obtained 
for the generic dimension a, and 73a (T (a)) is the Planck 
function. 

In the case of grains with dimensions larger than about 
~ 250A which remain cool because of their high thermal 
capacity and re-emit the energy absorbed from the incident 
radiation field into the FIR, a good approximation is to eval- 



uate the equilibrium temperature T (a) from the balance 
between absorption and emission 



J Qabs (a,X)J x dX = J Qats (a,X)B x (T(a))dX 



(8) 



where J\ is the incident radiation field. With the aid of T (a) 
we calculate the emission assuming that big grains of fixed 
composition and radius behave like a black body. Therefore, 
the emission coefficient is 



■ bin 



a Qabs (a, A) B x (T(a)) 



1 dn (a) 



mi 



da 



da (9) 



Emission from PAHs. Let us consider a photon of 
energy E = hv — hc/X impinging on a PAH molecule. The 
maximum temperature increase caused by the absorption of 
the photon is given by the solution of the equation 



he 

Y 



Cpah (T,a) dT 



(10) 



where Tmin ~ 5 K and Cpah (T, a) is the thermal capacity 
of a gr ain of dimension a at the temperature T. As pointed 
out bv IXu fc de Zotti £l989), T max (X , a) is essentially in- 
dependent from the choice of T m i n . In order to calculate 
Cpah (T, a), we use the nume rical fit of the iLeger et all 
l|l989h relationship proposed bv lSilva et al.l ll 19981) : 



T < 800 K 

800 K sC T «S 2100 K 

T > 2100 K 



Cpah (T, a) 

Cm ax 




(ll) 

where m = 9.25 • 10~ 4 n = 2 ■ 10" ' *, p = 0.58, and C ma x = 
3 [{No + N H ) - 2] k B (ILeger et al.lll989Tl . N c and N H are 
the number of hydrogen and carbon atoms in the molecule, 
and fes is the Boltzmann constant. 

Dealing with the number Nh of H atoms in a PAH 
an important issue has to be clarified. The PAHs in the 
ISM are subject to many physical and chemical processes, 
chief among which are the ionization and photo-dissociation 
caused by the absorption of a UV photon, electron recom- 
bination with PAH cations, and chemical reactions between 
PAHs and other major elemental species in the ISM. 

PAHs can be dehydrogenated (some of the C — H bonds 
can be broken), fully hydrogenated (one H atom for every 
available C site) or highly hydrogenated (some of the pe- 
ripheral C atoms bear two H atoms). The reason is that H 
atoms are less strongly bound (4.5eV) than carbon atoms 
(7.2eV) and therefo re can be m ore easily removed by ener- 
getic photons iLeger et al.lll989T) . 

Over the years many studies have been devoted to 
model the hydrogenation s ta te at various levels of sophisti- 
cation ]^^^^}^W96^ IVuong fc FoindbOOd) an d ie- 
centlv lLe Page et alJ (|2001l> and lLe Page et alJ l|2003h pre- 
sented a very detailed model for the charge and hydrogena- 
tion state of PAHs. The number of H atoms is given by the 
relation Nh = Xh ■ Ns where TVs is the number of avail- 
able sites for hydrogen to bind to carbon and \h is the 
hydrogenation factor. In a PAH with full hydrogen cover- 
age, the bond between TVs and TVh depends on nesting the 
hexagonal cycles into the molecule. These cycles can form (i) 
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an elongated low-density object, named catacondensed PAH 
(linear chains of n-cycles with general scheme C4 n +2^2n+4, 
like naphthalene with 2 cycles, anthracene with 3 cycles, 
and pentacene with 5 cycles); (ii) a compact object named 
pericondensed PAH described by the formula C§ p iH§ p , like 
coronene with p — 2, ovalene and circumcoronene or (iii) an 
untied object (biphenile). The more compact the PAH, the 
more stable it is. 

In agreement with this, to repre sent th e mai n 
characteristics of the in t erstel lar P AHs. | Omont] ||W86j), 
iDartois fc D'Hendecourll dl997fl and iLi fc Drainel (1200 ll) 
consider models of pericondensed PAHs, coronene-like, that 
are more stable against photodissociation by UV photons. 

From the formula for compact pericondensed PAHs 
with full hydrogen coverage, C 6p 2H(, p , we get Ns = 6p = 
(6 ■ 6p 2 )°' 5 = (6-/Vc) ' 5 - As the number Nc of carbon atoms 
is known, we can derive the number Ns of sites available 
for bonds. For example, for the coronene with p = 2 and 
Nc ~ 24, we have Ns = 12; if the coverage is full with 
Xh ~ 1, we have JV^r = 12. We are left with the hydro- 
genation \h factor to be determined. This will depend on 
the ISM conditions. For t y pical astronomical environments, 
ISieb enmorgen fc Krueeell il992T) have suggested dehydro- 
genated values from 0.01 to 0.5, wit hout excluding values 
as l ow as 0.01. More rece nt studies bv lLe Page et all (1200 ll) 
and lLe Page et all i2003fl adopt hydrogenation factors that 
depend on the PAH dimensions. Their analysis splits the 
population of PAHs in several groups. At the typical con- 
ditions of the diffuse ISM, their general conclusion is that 
very small PAHs (those with less than 30 C atoms) are ex- 
tremely dehydrogenated, larger PAHs with more than 30 
C atoms have full hydrogen coverage, whereas very large 
PAHs (with more than 150-200 C ato ms) can be high ly hy 



dro genated. As the model proposed bv lLe Page et al.l |2001 
and lLe Page et al] (|2003) is too much detailed for the pur- 
poses of_this_stu^ ! _th«_ff/C' ratios in usage here are taken 
from lLi fc Drainel (1200 ll) 



0.5 



H/C = \ 0.5/y/N c /25 
0.25 



N c «S 25 

25 < N c < 100 (12) 
N c 5* 100 



They yield N H = N s = (6Ac) ' 5 in the range 25 sC N c < 
100, which is typical of a pericondensed PAH with full hydro- 
gen coverage, while only for smaller PAHs there is a partial 
dehydrogenation. 

To calculate the emission of the PAHs we start defining 
the power irradiated by a molecule of dimension a at the 
temperature T 



S = f Tva 2 Q^ H (a, A) nB (A, T) dt — 
Jo 



T mi „ 



™ 2 Q P ab A s H (a, A) nB (A, T) ^ dT (15) 

F (T, a) 



where S = S{\ , A, a) and T max = T max (X , a). 

Under the effect of an incident radiation field I(X ) con- 
taining photons of any energy, the emission at the wave- 
length A per H atom becomes 



,,,, 1 f apAH dnpAH (a) 2r .PAH , ,s 
J (A) = / — —TraQ abs (a, A) x 



uh l„iom da 

PAH 

x r ^ x '™Q« bA " m x 

x f Tm ° X tvB (A, T) Cp /J' a } dT dX da (16) 



F(T, 



from which 



J (A) 



IX A 



a pAH dnpAM (a) 



da 



a2 Qabs H ( ci, A ) S I X , A, a ) dadX 



(17) 



The many types of PAHs can differ not only for the 
hydrogen coverage on peripheral carbon atoms, but also 
for the ionization state. Because of their low ionization po- 
tential (about 6 — 7 eV), on o ne hand the PAHs ca n be 
ionized by photo-electric effect ijBakes fc Tielenslll994h . on 
the other hand they can also acquire charge throug h colli- 
sions with ions and electrons iDraine fc Sutinlll987l) . Con- 
sequently, there is always so me probability that a PAH is 
in a charged state. Therefore. ILi fc Drainel {2001) have cal- 



culated and made available PAHs cross sections both for 
ionized and neutral PAHs. Since the ionization state will af- 
fect the intensity of AIBs, as we can see in Fig. in our 
calculations we have included the ionization of PAHs. 

To this aim, the charge-state distribution model of 
PAHs is taken from lWeingartner fc Dra inel (l200 1bl). who im- 
pro ved upon previou s attem pts Ty Tl3rarne^^Sutinl (Il987f) 
and iBakes fc Tielena dl994l) . Indicating with Z the charge 
state of a PAH molecule, in conditions of statistical equilib- 
rium we have the following relation 



F(T,a)= ™ 2 Qf£ H {a,X)nB(X,T)dX (13) 
Jo 

where B (A, T) is the Planck function. Once heated to a tem- 
perature T, the molecule cools down at the rate 



dt 



F(T,a) 



Cp{ T,a) ^ 

The energy emitted at the wavelength A by a PAH of dimen- 
sion a, as a consequence of absorbing a single photon with 
energy hc/X , is given by 



/ (Z) [J ve (Z) + ,h on (Z)] = / (Z + 1) J e (Z + 1) (18) 

where / (Z) is the probability for a PAH to possess the 
charge Z (either positive or negative), J pc is the photo- 
emission rate, Jj on is the positive ion accretion rate, 
and finally J e is the electron accretion rate. Following 
iDraine fc Sutml Jl987h . /(Z) is derived by recursively ap- 
plying eqn. 1181 to / (0). We obtain 



f(z)=f(o) n 



Jpe {Z 



1) + Jion (Z' - 1) 



Je (Z') 



(19) 
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Figure 1. Optical properties of neutral and ionized PAHs con- 
taining 40 atoms of carbon in the IR spectral region of the UIBs. 
The difference between neutral (dashed line) and ionized PAHs is 
soon evident: the ionization enhances the stretching modes C-C 
(6, 2 and 7.7 fim) and the bending mode C-H (8. 6 fim) , whereas i t 
weakens the stretching mode C-H at 3.3 fim iLi fc DraineH200ll) . 



for Z > 0, whereas for Z < we have 

Je (Z' + 1) 



/(Z)=/(0) J] 



(20) 



The constant / (0) is derived from the normalization condi- 
tion 



£/(*) 



(21) 



The expressions for J Q^^t^u>n 8.nd J e are taken from 
IWeineartner fc Drainel l|200HJ) . As both neutral and ion- 
ized PAHs are present, we indicate with \i — Xi ( a ) the 
fraction of ionized PAHs of dimension a, and with Qat^ H 
and Qabs AH the absorption coefficients of ionized and neu- 
tral PAHs. Including both types of PAHs into eqn. \17\ we 
obtain 



J (A) = 



+ 



nnhc 



-MA I A 



' h pAH dnpAH (a) 
aw da 



[oi 



I abs 



I Z AH (a,X)SioN(\,\,a) X r+ (22) 



a, A ) Sneu (A , A, a ) (1 — x0 dad\ 



3.2 The LID model 

We have also implemented in our code the physical and nu- 
merical treatment of t he inf rared emission of the ISM devel- 
oped bv lDraine fc"Tl l|200lf ) and lLi fc Drainel (|200J> . which 
is currently considered as the stat e-of-the art o f the s ub- 
ject and has already been used by ILi fc Drainel J2002bl) to 
study the IR emission of the SMC. Owing to its complex- 
ity, we will limit ourselves to mention here only a few basic 
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Figure 2. The energy distribution functions of silicate grains 
with dimensions of 10A, 22.5A, 35A, 56A, 89A, 14lA, 242A as in- 
dicated. All the distri butions have been c alculated u sing the ther- 
mal d iscrete model of lDraine &: Lil l200lt> . Following lLi &: Drainel 
i200 li) . we mark with a filled circle the first excited state. Po is 
the probability for the grains to be in the ground state that for 
small dimensions is high even if the distribution extends to high 
energy states. 



aspect s of the model. For all details the reader should re - 
fer to iDraine fc Lil feOOll) and ILi fc Drainel J200ll l2002bl) . 
In this model, made of three components, i.e. graphite, sili- 
cates and PAHs, the energy and temperature increase of the 
dust grains caused by the absorption of energetic photons is 
calculated by means of the density states of vibrational en- 
ergy levels. The task is accomplished in three possible ways: 
the exact statistical method, the thermal discrete approx- 
imation and the thermal continuous approach. The latter 
two make use of the same energy bins of the exact statisti- 
cal approximation, but differ in the way the transitions to 
higher (by absorption) and/or lower levels (by emission) are 
treated. Among the many improvements upon thermal mod- 
els made in the past, we call particular attention on (i) the 
detailed calculation of vibrational modes using the harmonic 
oscillator approximation, (ii) the much better estimate of the 
probability of grains to be in the fundamental state and (iii) 
the relationship between vibrational energy and "tempera- 
ture" of the grains. We have considered here both types of 
thermal model, leaving the exact statistical treatment aside 
because the enormous computational time required makes 
it difficult to implement it in spectrophotometric models of 
galaxies. In Fig. |5|we show the energy distribution functions 
of the silic ates for the th e rmal discrete model under the ef- 
fect of the iMathis et alJ il985t) radiation field. It is worth 
noticing how the probability for the small grains to be in 
the ground state is high, even if the distribution extends to 
high energy states. This accurate description of the ground 
state can not be obtained with the GDP model. 



© 2005 RAS, MNRAS 000, 000-000 



8 L. Piovan, R. Tantalo & C. Chiosi 



4 EXTINCTION AND EMISSION IN THE MW, 
LMC AND SMC 

As first step toward including dust emission and extinc- 
tion in the SED of young SSPs and galaxies with different 
shapes, different histories of star formation and chemical en- 
richment, different contents of gas, dust and metals at the 
present time, we seek to reproduce the average extinction 
curves of the MW, LMC and SMC and to match the emis- 
sion of their dusty diffuse ISM. As the MW, LMC and SMC 
constitute a sequence of decreasing mean metallicity, they 
provide the ideal workbench to calibrate the model for the 
ISM as a function of the metal content. Although includ- 
ing the effect of systematic differences in the metal content 
is of paramount importance, this aspect is no t alwa ys fully 
taken into account. For instance ISilva et al.l (1 19981) adopt 
for the diffuse ISM with different chemical compositions a 
model originally designed for the MW, and include the effect 
of different metallicities by simpl y scaling the abundanc e of 
dust with the metallicity (see also lDevriendt et alJll99fll for 
more details). 

The procedure we have adopted iteratively tries to si- 
multaneously satisfy the constraints set by emission and ex- 
tinction until self-consistency is achieved. It is also worth 
mentioning her e that our approach somehow improves upon 
the analysis bv lTakaei et alJ l|2003l ) . who did not test their 
dust model on the IR emission of the SMC. Unfortunately, 
we cannot test the model predictions for the emission of 
PAH s in the LMC because no data are available (see also 
iTakaei et afl l2003l. Nevertheless, making use of COBE and 
IRAS data, and coupling the fits of the emission and extinc- 
tion curves we can put some constraints on the flux intensity 
and the relative weight of the dust components. 

We present results obtained with the various models 
to our disposal, i.e. the GDP and LID model for emission 
and the WEI and our MRN-like distribution for the grain 
sizes. By comparing results from each possible combination 
we look for the one best suited to model the SED of SSPs 
and galaxies. 

4.1 The extinction curves of the MW, LMC and 
SMC 

The best fit of the extinction curves for the three galaxies 
under examination is derived from minimizing the y-square 
error fu nction for which we adopt the lWeingartner fc Draind 
fcOOlal) definition 

X 2 = E [A0bAV) 2 Am ° d(V)) (23) 

where A bs (Ai) /A b a (V) represent the observational data 
normalized to the V band, A mod (Ai) /A mod (V) are 
the values of the model, and o"i are the weights 
iWeingartner fc Draind l2001al) . Typically the number of 
test-points n is of the order of 100. We use as opti- 
mization techniq ue the Levenberg-Marquardt method (see 
IPress et al.l 11993 . for more details) combined with a sam- 
pling algorithm in bounded hypercubes and some manual 
pivoting of the parameters. In any case, it is worth notic- 
ing that the problem is very complicate because the multi- 
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Figure 3. Top Panel: The normalized distribution laws (num- 
ber density as a function of dimensions) for different types of 
grain. For the MRN model we show the graphite grains (thick 
continuous line), silicate grains (thick dashed line) and PAHs 
(thick dotted line). For the WEI model we display the carbona- 
ceous grains (thin dashed line) and silicate grains (thin continuous 
line). Bottom Panel: The extinction curve for the diffuse ISM 
of the MW obtained with a mixture made of graphite, silicates 
and PAHs. In this panel several cases are shown. First of all the 
total extinction for our MRN-like model (thick solid line) and for 
the WEI model (thin solid line) together with the observational 
data (same thin solid line). All the three lines are practically 
coincident. In addition to this for each case we show the con- 
tribution from the different types of grains: the silicates for the 
MRN and WEI models are indicated with the thick dashed and 
the thin dashed lines, respectively; the contribution of graphite 
in the MRN model and of carbonaceous grains in the WEI model 
are displayed by thick dotted and thin dotted lines, respectively; 
finally, for MRN model the contribution of PAHs is shown by the 
thick dot-dashed line. 



dimensional solution space is characterized by many local 
minima to which the solution may converge missing the true 
minimum. 

In Table 14.11 . we list the parameters we have found 
from the best-fit of the extinction curves of the MW, LMC, 
and SMC, respectively, using the MRN model. The WEI 
model is also shown for comparison. 

The MRN-like distribution given by eqn. JIJ has been 
used with n — 2 or n — 3 for graphite and silicates and n — 2 
for PAHs for the MW and LMC, whereas for the SMC a more 
complicated relationship has been adopted. Some parame- 
ters are kept fixed for all the models: i.e. the lower limits for 
the size distribution of PAHs and silicates which are set at 
the lowest values for which the cross sections are available, 
3.5A and 10A, respectively, and the upper limit of the sil- 
icates and graphite distribution fixed at 0.3 /jm and about 
1 fim. An important parameter is the upper limit of the PAH 
distribution, which is coincident with the lower limit of the 
graphite distribution: the two populations do not overlap 
and the PAHs are considered as the small-size tail of the car- 
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Figure 4. The same as in Fig.[3]but for the SMC. The data for 
the extinction curve of the SMC along the l i ne of s ight of the star 
AzV398 are taken from lGordon fc Clayton! <199d) . 
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Figure 5. The same as in Figs. |3] and ||] but for t he LMC. The 
data f or the extinction curve of the LMC are from lMisselt et alJ 
(1999). 



bonaceous grains. The transition size has been fixed at 20A 
thus allowing both PAHs and small, thermally fluctuating, 
graphite grains to concur to shape the SED in the MIR. The 
exact value of the transition size has in practice no influence 
on the extinction curve thanks to the similar UV-optical 
properties of PAHs and small graphite grains. The same con- 
siderations do not clearly apply to the IR emission, where a 
large population of PAHs wi ll produce stronger emission in 
the UIBs. On the contrary, in lWeingartner fc Drainel <l2001al) 



and lLi fc Drainel i200llh PAHs and graphite grains are in- 
cluded in the family of the so called "carbonaceous grains", 
where the smallest grains have the PAHs optical properties, 
the biggest grains have the graphite properties, and for the 
intermediate dimensions a smooth transition from PAHs to 
graphite properties is adopted. 

The final results for the extinction models, are shown in 
Figs.HH andglfbr the MW, SMC and LMC, respectively. 
The top panels show the scaled abundances of the various 
components in both the MRN- and WEI-models, i.e. PAHs, 
graphite and silicates in the former, carbonaceous and sili- 
cates in the latter. The bottom panels show the correspond- 
ing extinction curves, first separately for each component, 
then the total for both the MRN- and WEI-models, and fi- 
nally the observational data. The MRN- and WEI-models 
coincide and both are nearly indistinguishable from the ob- 
servational data. Furthermore, it is worth noticing that in 
order to obtain good fits for the MW, LMC and SMC the 
grain size distribution s need to be more co mplicate than the 
simple power law bv iMathis et al.1 lll977l) . Both in MRN- 
and WEI- models, the number of free parameters is about 
10. If a single power-law is adopted as in lMathis et alJ (Il977t) 
to reduce the number of parameters, we fail to get a good 
fit for all the three extinction curves. The worst situation is 
with the SMC. The absence of the bump in the extinction 
curve makes it difficult to get a good fit without increasing 
the number of adjustable parameters. The missing bump is 
indeed a strong constraint. Finally, our MRN-model best so- 
lution tends to yield abu ndances and overall results in c lose 
agreement with those bv lWeineartner fc Drainel {2001a). 



4.2 Emission by the diffuse ISM of MW, LMC 
and SMC 

Matched the extinction curves, we also need to test the over- 
all consistency of the model as far as the IR emission by dust 
and the comparison with the observational data for the SMC 
and MW are concerned . For LMC no data are available and 
the lLi fc Drainel <l200ll) model is taken as reference model. 

The two models of emission (GDP and LID) presented 
in Sects. |3~T1 and l3~2l must be coupled with the two models 
for the grain distribution (MRN and WEI) both of which 
fairly reproduce the extinction curves as amply discussed in 
Sect.Q 

The combination of WEI distribution w ith the LID 
model would be the id eal case to deal with (iLi fc Drainel 
l200ltlDraine fc Ll200ll) . However, we limit ourselves to con- 
sider only the cases GDP+MRN and GDP+WEI, because 
the GDP model is much faster than LID and therefore more 
suited to spectrophotometric synthesis. The LID model will 
be, however, considered as the reference case for comparison. 
The guide-line here is to cross check results for extinction 
and emission and by iterating the procedure to contrive the 
parameters at work so that not only unrealistic solutions 
are ruled out (that could originate by sole fit of the extinc- 
tion curve) but also additional information on the overall 
problem is acquired. 

The preliminary step to undertake is to ch oose the radi- 
ation field heating up the grains. We adopt the lMathis et all 
( 1983) interstellar radiation field (MMP) for the solar neigh- 
borhood of the MW, both for the whole MW and the LMC, 
whereas a slightly different radiation field is adopted for the 
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Table 1. Best fit parameters for the extinction curve of the diffuse ISM of the MW with Rv = 3.1, for the average extinction curve of 
the LMC and for the flat extinction curve of the SMC. For MW graphite and PAHs we adopt k = 2, whereas for silicates we use k = 3. 
For LMC silicates and PAHs we adopt k = 2, whereas for graphite we use k = 3. Finally, in the case of SMC, for PAHs we adopt k = 2, 
whereas for graphite we use k = 4 because we need a more complex set of power-law distributions. 
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Figure 6. Emission of the diffuse ISM predicted by models con- 
taining graphite, silicates and PAHs and comparison with the ob- 
servational data for the SMC Bar. The thick lines show the GDP 
emission model coupled with our MRN-like grain distribution for 
SMC. T he thin lines display th e thermal discrete LID emission 
model of lLi fc Drainel <200ll^ and lDraine fc LifcO Olft coupled wit h 
the WEI distribution for SMC iWeingartnerfc Drain j EoOl^ . 
The separate contributions from different grains are also shown: 
the solid lines are the total emission; the dotted lines are the emis- 
sion of silicates; the dot-dashed line are the ones of carbonaceous 
grains; finally the contribution of the stellar component is indi- 
cated by the dashed line. The in cident radiation field is the stan- 
dard MMP iLi fc Draindl2002rl . The data from COBE/DIRBE 
(diamonds) and IRAS (squares) have been kindl y provided b y Li 
(2004, private communication). See also lLi &: Drainel J2002b). 



SMC Bar. As this reg ion presents a large sp read i n radiation 
intens ities, following ILi fc Drainel <l2002bl) and iDale et ail 
i200ll) . the intensity of the radiation impinging on the dust 



grains has been described with a power-law distribution 
dN H /dU of MMP radiation fie l ds ILi fc Drainell2002 bl) . 

According to ILi fc Drainel £2002b) the intensity of the 
IR radiation emitted by the dust is 



- exp [-N^Vats (A)] 

N^a aba (A) 



d U^j x {U) (24) 



where, er a b s (A) is the total cross section of absorption, JVjf ' 
the total hydrogen column density and j\ (U) is the emission 
for the radiation intensity U, calculated as described in Sect. 

rm 

In Figs. HJ and we show the emission expected for 
the SMC, LMC and MW, respectively. Results for both 
GDP+MRN and LID+WEI are displayed (left panels) to- 
gether with the observational data, limited to MW and SMC 
(the Bar region); data for the LMC are not available. The 
comparison highlights the following: 

(i) In general the agreement between theory and ob- 
servational data is remarkably good. For the SMC, there 
is a small difference at about 10 /im: this is due to the 
differ ent distribution of silicates sizes, which according to 
ILi fc Drainel J2002b) may extend down to about 3.5A, in- 
cluding also nano-particles, for which, however, the cross 
sections are not yet available. 

(ii) No significant difference between GDP+MRN and 
LID+WEI can be noticed. 

(iii) In the case of the SMC the contribution of PAHs to 
the MIR is very weak: this follows from the absence of the 
bump in the average extinction curve that hints for a very 
low abundance of carbonaceous grains, while for the MW 
and LMC, the emission of PAHs becomes important. 

(iv) The only difference between GDP and LID to note 
is that the former overestimates the continuum underneath 
the 3.3 fim feature. This happens because GDP underesti- 
mates the probability of grains being in the ground state. 
As a consequence of it, the grains are generally hotter and 
thus emit more energy at shorter wavelengths. 

(v) To get a deeper insight on the performance of the 
GDP model, we consider the case GDP+WEI and apply it 
to predict the emission from MW and LMC. In such a case, 
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Figure 7. Left Panel: Emission of the diffuse ISM predicted by models containing graphite, silicates and PAHs appropriate to the 
case of the LMC. The thick lines correspond to the GDP emission model with our prescriptio n for MRN-like d i stribu tion law suited 
to the LMC. Th e thin lines are the thermal di screte LID emission model of lLi fe Drainei 1200 J) and lDraine fe lH r 
distribution law IWeingartner 



1200 ID with the WEI 

The separate contributions from different grains are also shown: the solid lines are the 



DrainelEooTal 

total emission; the dotted lines are the emission of silicates; the dot-dashed lines are the ones of graphite grains; finally, the short dashed 
lines are the total emission of carbonaceous grains, graphite plus PAHs. The incident radiation field is the standard MMP. Right Panel: 
the same as in the Left Panel but in all cases the distribution law of [Wcingartncr & Drain^ l2001aTl is adopted. 
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Figure 8. Left Panel: Emission of the diffuse ISM in the MW towards high Galactic latitudes (|6| 25°). The emission model takes 
into account graphite, silicates and PAHs. The thick lines correspond to the GDP emission mode l with our prescrip tion for MRN-like 
distribution law suited to the MW. The thin lines arc the thermal discrete LID emission model of lLi fc Drainei 1200 ll) and lDraine fc Lil 
(2001) with the WEI distribution law IWeingartner fc Dr ainc 2001aj). The separate contributions from different grains are also shown: 
the solid lines are the total emission; the dotted lines are the emission of silicates; the dot-dashed line are the ones of graphite grains; 
finally, the short dashed lines are the total emission of carbonaceous grains, graphite plus PAHs. The incident radiation field is the 
standard MMP. The observational data from DIRBE (diamonds) and FIRAS (squares) have been kindly provided by Li (2004, private 
communication'! . See | L i fc Dr ainei 1200 ll) . Right Panel: the same as in the Left Panels, but in all cases the same distribution law of 
IWeingartner fc DrainaT2001al) has been adopted. 
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both GDP and LID models use the same grain distribution 
so that the sole effects of the emission models can be singled 
out. The results are shown in the right panels of Figs. [7| 
(LMC) and 03 (MW). Once more, not only the GDP model 
yields results that agree with the observational data, but 
also with LID. There are two marginal differences: one at the 
3.3 (j,m feature and the other at about 50/im where the flux 
is slightly underestimated. The reason is always the ground 
state probability: in GDP grains tend to be hotter and to 
emit more at short IR wavelengths thus slightly shifting the 
emission from 50 /im to shorter wavelengths. Finally, using 
in both GDP and LID models the WEI distribution, the 
small differences noticed for LMC at 100/im also disappear. 

Basing on this combined analysis of emission and ex- 
tinction, we decided to adopt GDP+WEI to describe the 
dusty ISM and to model the SED of young dusty SSPs and 
of galaxy spectrophotometr ic models tha t will be presented 
in the companion paper bv lPiovan et alJ J2005f) . 

There is a final consideration about our choice for the 
size distribution law. Even if both WEI and MRN lead to 
good fits, it is only the emission model that determines the 
correct solution. It is worth recalling that it is not possi- 
ble to constrain the C abundance only by fitting the ex- 
tinction curve which simply provides an upper limit (see 
IWeingartner fe Draind|20 01a. for more details). With the 
MRN model it is much more difficult to obtain good fits of 
the extinction curves at varying the C abundance and to 
couple extinction and emission. In contrast, with the WEI 
model, which considers be as a free parameter, the goal can 
be easily achieved. Therefore, the WEI distribution law is 
perhaps best suited to model the SEDs of young dusty SSPs 
and galaxies. 



5 POPULATION SYNTHESIS FOR SSPS 

The monochromatic flux of a SSP of age t and metallicity 
Z at the wavelength A is defined as 



sspx it, Z) 



Mu(t) 



fx(M,t,Z)$(M)dM (25) 



where fx (M, t, Z) is the monochromatic flux emitted by a 
star of mass M, age t, and metallicity Z; <E> (M) is the IMF; 
Ml is the mass of the lowest mass star in the SSP, whereas 
Mu(t) is mass of the highest mass star still alive in the 
SSP of age t. For the IMF we adopt the ISalpeterl lh955Tl 
law expressed as dN/dM = $(M) = AM~ X where £=2.35 
and A is a normalization constant to be fixed by a suitable 
condition (SSPs for other choices of the IMF can be easily 
calculated) . 

In our study we adopt the isochrones bv lTantalo et all 
l|l998fl (anticipate d in the data base for galaxy evo- 
lution models by iLeitherer et al The underly- 
ing stellar models are those of the Padova Library (see 
iBertelli et all Il994 for more details). The initial masses 
of the stellar models go from 0.15 to 120 M©. The fol- 
lowing initial chemical compositions have been consid- 
ered: [Y=0.230, Z=0.0004], [Y=0.240, Z=0.004], [Y=0.250, 
Z=0.008], [Y=0.280, Z=0.02], [Y=0.352, Z=0.05], and 
[Y=0.475, Z=0.1], where Y and Z are the helium and metal 
content (by mass). 



The s tellar spectra i n usa ge here are taken from the 
library of iLeieune et alJ <ll998l) . which stands on the Ku- 
rucz (1995) release of theoretical spectra, however with sev- 
eral important impleme ntations. For T e ff < 3500 K the 
spectra of dwarf stars bv lAllard fc Hausclnkhj j[l995j) are in - 
clud ed and for gia n t star s the sp ectra by lFlTfksetHfll dl994h 
and iBessell et al] Jl989l 1991) are considered. Following 
iBressan et alJ Jl994h . for T eff > 50000 K, the library has 
been extended using black body spectra. 

The SSP spectra have been calculated for the same 
ranges of ages and metallicities of the isochrones. 



6 SPATIAL DISTRIBUTION OF YOUNG 
STARS, GAS AND DUST 

The first step to undertake is to specify the relative dis- 
tribution of young stars and dust. From the observational 
maps of MCs it is soon evident that the situation is very 
complicate: there are many point sources of radiation (the 
stars) enshrouded by dust, that are randomly distributed 
across clouds of irregular shape. Inside a cloud the tem- 
perature varies locally depending not only on the distance 
r from the cloud center but also on the distance from the 
nearest star. In other words a dust grain not particularly 
close to a star will have a temperature T (r) determined by 
the local, average interstellar radiation field (ISRF) , whereas 
a dust grain at the same distance r from the center, but 
close to a hot star, will exper ience a hotter temperature 
llKrugel fc Siebenmorgenl Il994f) . All this implies that the 
spherical symmetry is broken, thus drastically increasing 
the complexity of the whole problem. The assumption of 
spherical symmetry is the only way to handle this kind of 
problem at a reasonable cost. However, even with spherical 
symmetry, the spatial distribution of stars, gas and dust can 
be realized in many different ways. Finally, there are differ- 
ent techniques to deal with the radiative transfer across the 
den se medium of the MCs. 

iKrueel fc Siebenmorgenl lll994fl analyzed in a great de- 
tail three possible spherically symmetric configurations: first 
the case in which hot stars enshrouded b y dust (the so-called 
"hot- spots") are randomly distributed llKrugel fc Tutukovl 
ll978Tl : secondly the case of a central point source 
(the "point-source" model) (see Rowan- Robinson! Il98d : 
iRowan-Robinson fc Crawford! Il98£ ) and finally the case 
with an extended central source Jsiebenmorgenlll99 if) . 

The hot-spots description (or anyone of the same type) 
would be closer to reality and the one providing more physi- 
cally sounded results. The back side of the coin is that com- 
pared to the central source case, the problem gets soon com- 
plicate from the point of view of numerical computations. 
For this reason we take an hybrid case in between the hot- 
spots and central point-source descriptions. We simulate a 
young stellar population embedded in a spherical cloud of 
gas and dust by assuming that stars, gas and dust follow 
three different King's laws each of which characterized by 
its own parameters: 



Pi (r) = po,., 



1+1 — 



(26) 



where the index i can be s (stars), d (dust) and g (gas) 
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Figure 9. Spectra of dusty SSPs of metallicity Z = 0.02 calcu- 
lated with DUSTY (central point source - dotted lines) and with 
the ray tracing method (distributed source - continuous lines). 
They are represented at three different ages, from top to bottom: 
10 Myr, 100 Myr and 1 Gyr. 

rc,s,r c , a and r c ,d are the core radii, whereas 7 9 ,7d and 7 a 
are the exponents of the distributions of gas, dust and stars. 
As dust and gas are mixed together, they are described by 
the same parameters, r Cj£ j = r Cj9 . The exponents are simpl y 
chosen to be 7 S = ja = 7a = 1.5. As in Tak agi et, alJ |2003), 
we introduce the parameter r\ = r c ,d/r c , s , equal t o the r atio 
between the two scales. Following ICombes et alJ a!995f) we 
adopt the relation log (rt, s /r c ,s) = 2.2 between the tidal ra- 
dius rt, s and the scale radius of stars r c , s and simply assume 
rt,s = r t , g = r t ,d = r t . Denoting with Rmc the radius of a 
generic MC, and imposing rt = Rmc , there will be no d ust 
and gas for r > Rmc, i.e. for r > rt IITakagi et aij 2003). 

To fully understand the differences between a central 
point-source and a spatially distributed source it may be 
worth of interest to compare SSPs calculated with the two 
schemes. The SSPs for the spatially distributed source are 
the ones presented here. They are calculated using the ray 
tracing technique to be explained in Sect. 17.11 For SSPs 
based on the central point-source description we adopt 
the code DUSTY 1 to solve the radiative transfer equation 
Jlvezic fc Elftzurlll997f) . Even if DUSTY is not able to han- 
dle thermally fluctuating grains and PAH emission, it is fully 
adequate to show the difference arising from different geome- 
tries. 

In DUSTY only two quantities must be specified to ob- 
tain a complete solution: the total optical depth t\ at a 
suitable reference wavelength, and the condensation temper- 
ature of the dust (T s ) at the inner edge of the shell. All other 
parameters are defined as dimension-less and/or normalized 
quantities. We h ave fixed T s to a suitable value, following 
the suggestion bv lSilva et al.l |l998) who made use of central 
point source approach to model MCs. 

In Fig.Hwe compare young dusty SSPs calculated both 
with DUSTY (central point-source) and the ray tracing 

1 The code is kindly made available by the authors at 
http://www.pa.uky.edu/~moshe/dusty. 



technique (distributed source). Both SSPs have the solar 
metallicity, the same optical depth and matter profile of the 
cloud, and the same dust properties. The geometry is the 
origin of the lack of stellar emission for the dusty SSPs cal- 
culated with the central point-source model. All the stars 
are indeed embedded in a medium of high optical depth so 
that no stellar light can escape from it. Even if in some cases 
of highly embedded sources this could be a good approxi- 
mation, in general it is not realistic, because we can easily 
conceive that in a real environment with ongoing star for- 
mation, young stars can be formed everywhere even near 
the edges of the MCs. Such stars are therefore less obscured 
by dust. The ray-tracing method with the stars distributed 
into the cloud can easily handle this situation thus allowing 
a fraction of the stellar light to be always detectable. This 
is a net advantage with respect to the central point-source 
approximation. As a matter of facts, there is observational 
evidence for contributions to the UV-NIR fl ux from bursts 
of you ng star s as discussed bv lTakagi et alJ (see the discus- 
sion by 20Q3). The central point-source model is not able to 
simulate properly this contribution, which is instead nicely 
described by young stars being distributed all across the 
cloud. Anyway, as the point-source model can reasonably 
simulate highly obscured and concentrated populations of 
stars, we plan to introduce also this case in our library of 
SEDs. 

7 THE RADIATIVE TRANSFER PROBLEM 
7.1 The ray tracing method 

To solve the radiative transfer equation we use the robust 
and simple te chnique otherwise known as the "ray trac - 
ing" method feand fc Grindlavl Il98sl: iTakagi et alll2003l) . 
The method solves the radiative transfer equation along a 
set of rays traced throughout the inhomogeneous spherically 
symmetric source. Thanks to the spherical symmetry of the 
problem, it is possible to calculate the specific intensity of 
the radiation field at a given distance from the center of 
the MC by suitably averagin g the intensities of all th e rays 
passing through that point feand fc Grindlavl 11985). The 
geometry of the problem is best illustrated in Fig. HOI which 
shows a section of a spherically symmetric MC. We con- 
sider a discrete set of N concentric spheres (circles on the 
projection) and an equal number of rays indicated by the 
parameter j running from j=l to j=N. The specific inten- 
sity is calculated at all intersections of the generic ray j with 
the generic circle (sphere). The number of intersections i in- 
creases from 1 for j=N to 2(N — 1) + 1 according to the 
rule i=2(N — j) + 1. The equation of radiative transfer for 
specific intensity I\ (y, x) along a ray can be written as 

dh ^.' X) = -«fl (r) °e*t (A) [h (y, x) - S X (r)] (27) 

where y is the impact parameter, i.e. the distance of the 
generic ray j from the ray passing through the center, and 
x is the coordinate along the ray; uh (r) is the number of 
hydrogen atoms per cubic centimeter; a ex t (A) is the cross 
section of extinction for hydrogen atoms defined in eqn. JSJ; 
finally, S\ (r) is the source function at a given radius. We in- 
clude the metallicity dependence for the composition of the 
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Figure 10. Two dimensional section of a dusty MC. A set of TV 
rays is traced each one tangent to one of N concentric spherical 
surfaces of a radial grid with TV = 10 and thicker spacing toward 
the center of the cloud. The specific intensity is calculated for each 
ray at each point when the ray penetrates one of the surfaces. 
The j-th ray cross TV — j surfaces in 2 (TV — j) + 1 points, j/j is 
the impact parameter of the j'-th ray. 

MCs (on the notion that a generation of newly born stars 
shares the same metallicity of the surrounding MCs) as fol- 
lows: for metallicities Z ^ 0.015 we use the extinction curve 
of the MW for dense regions, characterized by a reddening 
parameter as high as is Rv = 5.5; for metallicities in the 
interval 0.005 ^ Z ^ 0.015 we adopt the extinction curve of 
the LMC; finally, for metallicities lower than Z ^ 0.005 we 
use the typical extinction curve of the SMC. 

In order to solve eqn. (1271 . following iRowan-Robinsonl 
(1980) we split the specific intensity in three components: 
I x = + l[ 2) + if ) , where 5 is the light from the stellar 
source, J^ 2 ' is the light scattered by the grains, and, I x 3 ^ 
is the radiation coming from the thermal emission of the 
grains. It is possible to calculate the contribution from each 
component, if the optical properties of the dust grains do 
not depend on their temperature. T his is indeed the cas e 
of the optical coeffici ents taken fromlDraine fc Lee) lll984) . 
lLaor fc Prairie! Jl995t and lLi fc Draind d200ll) . As we are in 
spherical geometry only a radial grid needs to be specified. 
For each impact parameter yj = rj, I^ 1 ' is calculated and 
stored at each intersection of the j-th ray with the spheres 
(circles) (as shown in Fig. llOII . where the stellar source func- 
tion (r) is taken from lTakagi et all i200ot) . 

Once the matrix Iy' (yj,x i+ i) is calculated, it is pos- 
sible to obtain the mean intensity J^' (rj) at a given ra- 
dius performing an integration over the cosine grid \ij as 
displayed in Fig. 1111 We have now to calculate the contri- 
bution of l[ 2) and lf> to J A . The source function of the 
n-times scattered light is 

(^) = ^ / <M™> l( xi-i (r, m) dQ (28) 

(2) 

where I\ n _ 1 is the associated specific intensity. The func- 
tion (f>(jn,Ti) is the phase function, which yields the dis- 
tribution of photons scattered in different directions. For 
the sake of simplicity we assume isotropic scattering so that 




Figure 11. Two dimensional section of a dusty MC. A set of 
TV spherical surfaces of a radial grid with TV = 10 is shown. The 
average intensity is calculated for the point j' by using all the rays 
from 1 to j' — 1, each rays weighted for the cosine of the angle 
with the radial vector. For each ray we have two contributes, one 
for each direction, Xj_ and Xj, . 



<f>(m, n) = 1. The total intensity of scattered light is ob- 
tained by adding all scattering terms I x 2 ^ = Yln^Xn- We 
start from I x \ which can be considered as the 0-times scat- 
tered light, to obtain the source function of 1-times scattered 
light (Ta kagi et alJl2005T) and the average specific intensity 

(2) 

for each spherical concentric shell. We derive S x n (r) and re- 
peat these steps, calculating the scattering terms of higher 
order. By taki ng a sufficiently l arge number of steps, six 
are sufficient dTakagi et al.1 [20031. the energy conservation 
is secured. 

Finally, we need the source function for the infrared 
emission of dust. We take also dust self-absorption into ac- 
count: for high optical depths such as the ones of MCs, the 
infrared photons can be absorbed and self-absorption has to 
be considered. An iterative procedure is required until the 
energy conservation is reached: using the mean intensities 
of direct and scattered light, (r) and (r), the dust 
source function for the 1-th iteration (r) together with 
the corresponding average intensity [ (r) are obtained. 
Starting from (r), we iterate the procedure to get the 
thermal emission of dust J^ n (r) by means of the the local 
radiation field jf 5 (r) + (r) + J^_j (r). The dust emis- 
sion is calculated as described in Sect. 13.11 with the GDP 
model. Therefore, the source function of dust is given by 

S$ (r) = (jr a " + jl' 3 + f, AH ) (29) 

where j s x mal1 , j" s and jf AH are the contributions to the 
emission by small grains, big grains and PAHs, respectively. 
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7.2 Optical depth of the cloud 

Key parameter of the radiative transfer is the optical depth 
t\ defined as 



t\ = 



kxpds 



(30) 



where k\ is the extinction coefficient, p is the density of the 
matter and S is the thickness of the cylinder of matter along 
which we integrate. Eqn. I30H can be recast as follows 



the density law for a MC given by eqn. (12611 over its volume 
and using the identity 

r-RMc/ r c, g x 2 
Mmc = 4tt po, g (r c , g ) (x + a . 2 )-i.5 dx ( 35 ) 

with x — r/r c ,g- 

The above expression for the optical depth depends only 
on suitable constants and the dimension-less integral which 
is a function of rj. Denoting with 



t\ = 



n H (r) ax (r) dr 



(31) 



where nn (r) is the number density of H atoms and o\ (r) 
is the extinction cross section (scattering plus absorption). 

In order to calculate the thickness of the shell, we need 
an est imate of the mas s Mmc and the radius Rmc of the 
cloud. Isilva et al1<ll99 8) adopt a "typical" MC of given mass 
and radius. However, there is plenty of evidence that MCs 
span large ranges of masses and radii in which the low mass 
and small radius ones are by far more numerous. Therefore, 
in a realistic picture of the effects of dust clouds, the mass 
and radius distributions of these latter should be taken into 
account thus requiring that large ranges of optical depths 
are considered. We start calculating the optical depth for 
a cloud with a generic mass Mmc and radius Rmc- From 
the identity p g (r) = p H (r) + pHe (r) + pz (r) with p z (r), 
Ph (r), pHe (r) densities of metals, hydrogen, and helium, re- 
spectively, after putting pn (r) into evidence and neglecting 
pz (r) Iph (r) we get: 



Ph (r) 



Pa ( r ) 



PHe (r) 
Ph (r) 



Po,g- 



-~ig 



PHe (r) 

Ph (r) 



(32) 



We can now reasonably assume that the hydrogen and 
helium distributions have the same scale radius, so that 
pHe (r) I pn if) = pHe.o/ pH,o does not depend on the co- 
ordinate r. Dividing eqn. (1321 by the mass of the H atom 
run and supposing that o\ (r) does not depend on the co- 
ordinate r along the cylinder of matter we derive 



-, i pHe,0\ ,, 



1 + 



-7s 



dr (33) 



Setting z = r/r c , 9 and introducing the parameter r\ — 
r c ,d/r c ,8 = r Ctg /r c . a - with the ratio of the cloud radius 
Rmc to the scale radius of the gas component given by 
RMc/r c , g = 10 2 ' 2 • (r c ,s/r c ,s) = 10 2,2 /r/ - and 7 9 = 1.5 we 
obtain 



uxPo.qRmc 



10 2 - : 



(1 + z 2 ) 1 ' 6 dz (34) 



The normalization constant po, g is derived by integrating 



*(»?) = 



{i+z 2 y x - h dz 



10 2 jo 

and K = 1/ [(1 + pHe.o/pn.o) m H ], we finally get 
t\ = KcfxPq^Rmc $ {v) 



(36) 



(37) 



The parameter rj describes the dust distribution: for 77 = 
1000 we have a nearly uniform distribution of dust in which 
stars are embedded, whereas for 77 = 1 stars and dust in 
the MCs have the same distribution. Typical values of "I> (rj) 
are $ (77) = 0.0063, 0.063, 0.534, 0.988 as 97 goes from 1 to 
10, 100, 1000, respectively. Following iTakari et afl (|2003h 
who reproduced the central star forming regions of M82 
and Arp220 using a nearly uniform distribution of dust 
(rj = 1000), we adopt the same value for this parameter. 
Looking at the optical depth of the cloud we have just de- 
fined, there are other considerations to be made about the 
geometry of the system. Both the mass Mmc, via the central 
density po, a , and the radius Rmc of the cloud are present 
in eqn. 13711 . The first problem to deal with is that masses 
and radii of real MCs span large ranges of values. This prob- 
lem will be addressed in detail in Sect. 17.31 Second, in all 
the relationships we have been using in so far to define the 
final expression for the optical depth, the total mass of the 
stars embedded into the MC does not play any role. This 
is obvious because the optical depth depends only on the 
dust-fog in which the stars-candles are embedded, and not 
on the candles themselves. So, which star mass M* has to 
be considered for a cloud with mass Mmc and radius Rmc7 
There is no easy answer to this problem. We start consider- 
ing that our MCs spectra will be eventually used to build in- 
tegrated spectra of galaxies. Therefore, we need normalized 
MCs spectra to be convolved with the SFH of the galaxy. 
The simplest way we can achieve the goal is to consider as 
candle a SSP of M* = IMq. We keep fixed the optical depth 
of the cloud calculated with eqn. 13 7H using the real dimen- 
sions of the cloud and then w e re-scale the radiu s of the 
cloud according to the relation jTakaei et alJl2003|) 



Rm 



WRs 



M k 



10 U M G 



(38) 



where R' M c ( m kpc) is the normalized radius corresponding 
to the adopted mass M t (for an arbitrary star mass M*, 
the radius Rmc ha s to be suitably scaled). As noticed by 
Taka gi et alJ <|2003T) . when the optical depth and the scaling 
parameter R a are kept fixed, the relative shape of the SED 
remains unchanged for different masses M*. The only effect 
of varying is to produce a higher energy output, but the 
SED remains the same. 
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7.3 The mass spectrum of MCs 

As already mentioned, there is observational evidence that 
the MCs span large ranges of masses and radii. Let us as- 
sume now that masses and radii of molecular clouds obey the 
distribution laws dN/dMiac = AmMm'c and dN/dRMC = 
AnR^f c where the normalization constants An and Am can 
be determined integrating over the distributions once the 
lower and upper limits are specified. We split the mass and 
radius intervals in a discrete number of bins, and indicate 
with n and m the number of the mass and radius bins, re- 
spectively. The normalized number AT; (Mmc) /N of clouds 
into the i mass bin of n bins and the number of clouds 
Ni (Raw) /N into the i radius bin of m bins can be easily 
obtained. They are independent of the normalization con- 
stant An and Am- The sum over all the bins is of course 
normalized to unity. The fraction of molecular clouds with 
mass belonging to the i-th bin of n bins and radius belonging 
to the J-th bin of m bins will be 

/(M^) = ^^ (39) 

According to eqn. I37H . for r\ — 1000 the function $ (77) ~ 1 
and the central density po, g is nearly equal to the mean 
density because dust is almost uniformly distributed, i.e. 
Po,g ~ 7> g - Recalling that p g = 3Mmc/ (4-kR mc ), we ob- 
tain t\ cc Mmc I R\ic- This means that for each ratio 
Mmc/R\ic a different value of the optical depth is ex- 
pected. Furthermore, splitting the range of optical depths in 
k bins, it may happen that different combinations of Mmc 
and Rmc yield a similar value of t\. As a consequence of 
this, different pairs of bins may give their fractional 

contribution to the same bin r of the optical depth. 

The limits of the mass and radius distri butions can be 
inferre d from the observational surveys of ISolo mon et alJ 
J1987D of galact ic M Cs an d the analysis made by 
lElmegreen fc Falearond <ll996T) . Typical values are: R up ~ 
70 pc, Ri ow ~ 5pc, M np ~ 3x 10 6 M© and finally M iow ~ 4x 
10 4 Mq. The exponents of the mass and radius distribution, 
ol and P respectively, are taken from lElmegreen fc Falearonel 
( 1996). We have /3 ~ 3.3 ± 0.3 for the radii distribution and 
a — — 1.8 for the mass distribution. These power-laws imply 
that not all the optical depths r(i,j) will have the same 
weight, but those falling in suitable fc-bins will weigh more 
so that the whole range can be reduced to these values. 
This is clearly shown in Fig. I12l where the discrete function 
f(Mi,Rj) is plotted over 15 radius bins and 15 mass bins 
and the /(l, 1) bin weights much more than the other ones. 
However, there is a strong dependence on the lower limits of 
the distribution laws we have adopted, i.e. the one holding 
for the MW, and it is likely that in different environments, 
such as strong star forming galaxies, the distribution laws 
may be different. Therefore, a realistic model of galaxies 
should include MCs with an ample range of optical depths 
and also different laws for the mass and size of the MCs. For 
the purposes of this study, we limit ourselves to the distribu- 
tion laws holding for the MW. With these distribution laws 
and using 15 bins, the /(l, 1) bin weighs about 50% of the 
total, as shown in Fig. 1121 Thanks to this, for the time being 
we can simplify the problem and take only one optical depth 
as a free parameter representative of the whole range, both 




j-th Radius bin i-th Mass bin 

Figure 12. Three dimensional plot of the function / (Mi, Rj) for 
15 mass bins and 15 radius bins. The first bin corresponds to the 
lowest mass and radius of the distributions. 

in sta r forming regions (Sect.|3J and galaxies dPiovan et alJ 
2005). 



8 SEDS OF YOUNG SSPS SURROUNDED BY 
MCS 

In this section, we present the SEDs of young SSPs in which 
the effects of dust around are included according to the pre- 
scriptions of Sect. |7| and compare them with the standard 
case in which dust is neglected. For each metallicity the anal- 
ysis is limited to the most probable bin /(l, 1) of the mass 
and radius distribution of MCs. The optical depth of this 
bin is about tv = 35 for the solar metallicity, tv = 10 
and tv = 5 for the LMC and SMC metallicities (extinc- 
tion curves). The scaling parameter for the MC dimensions 
is R = 1, while the total carbon abundance per H nucleus 
in the log-normal populations of very small carbonaceous 
grains, b c , has been fixed to a value slightly above the mean 
value of the a llowed range. Finally, the trea tment of ion- 
ization is as in lWeingartner fc Drainel <l2001blh More details 
will be given in Sect. 18.11 describing the library of young 
SSPs with dust we have calculated. 

We begin showing in Fig.EUla SSP with Z = 0.02 (solar 
metallicity) and age of 30 Myr, both in presence and absence 
of dust. The main effects of dust are soon clear: the cloud in 
which the young population is embedded absorbs the stellar 
radiation and returns it in the IR. We may also note here 
the different contributions by the various components of the 
specific intensity: direct light, scattered light and dust emis- 
sion. Passing to discuss the effect of the metallicity, for the 
sake of brevity we will show results only for Z = 0.02 (Fig. 
ITHl. Z = 0.008 (Fig. ESI, and Z = 0.004 (Fig. [Toll. 

As expected, the SED shift toward the IR is stronger 
for solar or LMC metallicities, because of the higher optical 
depth of the MCs (all the other geometrical parameters of 
the cloud being kept constant). For a given metallicity the 
peak in IR emission of cold dust shifts toward shorter wave- 
lengths at decreasing age of the SSP, because the intensity 
of the radiation field gets higher. For each shell of the radial 
grid in which a MC has b een split, we calculate t he ion- 
ization state of PAHs. The IWeingartner fc Drainel fcOOlbl) 
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Figure 13. SED of a SSP of 30 Myr and solar mctallicity with 
and without dust. In the plot we show the original spectrum of the 
SSP (thin solid line), the total dusty spectrum (thick solid line), 
the contribution of direct light (dot-dashed line), the contribution 
of scattered light (dashed line) and, finally, the emission of dust 
(dotted line). 




Figure 14. SEDs of SSPs of mctallicity Z = 0.02 without dust 
(dotted lines) and with dust (solid lines). The spectra are obtained 
solving the radiative transfer problem by means of the ray tracing 
method. The SSPs are shown at four different ages, from top to 
bottom: 3 Myr, 10 Myr, 100 Myr and, finally, 1 Gyr. 



model yields SEDs of very young SSPs in which the PAHs 
are almost fully ionized and the profile of aromatic infrared 
bands of PAHs (AIBs) is very close to that for the ionized 
cross sections, as shown by Figs. 1141 1151 and 1161 At increas- 
ing age and hence decreasing radiation field, the profile of 
AIBs shifts to that of neutral cross section because more 
and more PAHs become neutral. For the lowest metallici- 
ties the contribution of PAHs becomes small because of the 
small percentage of carbonaceous grains in the adopted ex- 




Figure 15. The same as Fig. ITlbut for the metallicity Z = 0.008. 




Figure 16. The same as Fig.Elbut for the metallicity Z = 0.004 
and the youngest age that is 3.5 Myr. 



tinction curve (the one for SMC). Furthermore, at decreas- 
ing age, AIBs become weaker due to the combined effects 
of the lower optical depth and higher radiation field which 
cause emission by hot dust in the same wavelength interval 
of AIBs. 



8.1 A library of young SSPs surrounded by MCs 

Dealing with theoretical spectra of young SSPs embedded 
into dusty MCs is a cumbersome affair as a large number of 
parameters is involved. First, we have those related to dust 
and its composition (still uncertain), together with those of 
the optical properties and distribution laws of the grains. 
Second, we have the geometrical parameters introduced to 
solve the radiative transfer in a medium with sources dis- 
tributed in it. Although one could simply take a sort of pro- 
totype MC to be inserted in population synthesis studies, 
it is more physically grounded and close to reality to create 
and explore in some detail the space of parameters. 

First of all, let us discuss the optical depth t\ and the 
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scaling factor R of MCs (the latter fixing their s i ze). F or 
more details about R and t\ see also lTakagi et all J2003) . 

T\: the optical depth, see eqn. 1371 . is the main pa- 
rameter of the radiative transfer problem. Its effect on the 
spectrum emitted by the MC is simple to describe: if the 
optical depth is high, more energy is shifted toward longer 
wavelengths. This amount of energy first quickly increases 
at increasing optical depth, then becomes less sensitive to t\ 
and even tends to flatten out for high t\ . This can be seen in 
Figs. I14) and ll5l where even if the optical depth is different, 
yet it is high enough to show the above effect. The opti- 
cal depth depends on the mass and radius of the MC (the 
geometrical parameters) and different distribution laws for 
masses and radii lead to different optical depths. The ideal 
case would be to set up a library of SSPs covering an ample 
range of optical depths. Because of the high computational 
time required, the present release of the library is limited to 
two values of the optical depth for each metallicity, namely 
tv = 35 and ry = 5. 

R: the scaling parameter R has been introduced because 
the mass M* contained inside a MC has been normalized to 
1 Mq, in view of the use of these SSPs in evolutionary popu- 
lation synthesis models for galaxies. R links the mass of the 
sources of radiation to the dimension of the cloud, as de- 
scribed in Sect. 17.21 There is no effect of varying R over the 
part of the SED given from the absorbed stellar emission. 
The only effect of R is to change the position of the FIR 
peak due to dust emission. A ideal MC scaled with a larger 
value of R will have a lower temperature profile of the grains 
because of the bigger dimensions. Accordingly, the FIR peak 
of dust emission will simply shift to longer wavelengths. This 
effect is entangled with the age effect, because for older pop- 
ulations the radiation field will be weaker and the FIR peak 
will undergo a similar shift to longer wavelengths. 

Let us now examine the parameters controlling the MIR 
region of the spectrum, i.e. the intensity of PAH emission 
and the AIB profiles. 

be- this parameter is the total C abundance (per H 
nucleus) in the two log-normal distribution laws of very 
small carbonaceous grains. Even if other distribution laws 
for the grain dimensions can be us ed, our preference goes to 
the IWeingartner fc Drainel i2001ah law because it has been 
proved to provide excellent fits to the extinction curves and 
it contains the contribution of very small grains to carbon 
abundance as a free parameter. Only an upper limit to the 
parameter be can be derived and all the values lower than 
or about equal to this upper limit are possible (see Sect.^. 
However the exact value of be bears very much the PAH 
emission in the MIR. At increasing be, the AIBs are more 
pronounced with higher flux levels, while the flux in the UV- 
optical region and the shape of the absorbed stellar emis- 
sion remain unchanged as the extinction curve is exactly the 
same for different values of be- At decreasing be, the emis- 
sion in the FIR slightly increases (the total energy budget 
has clearly to be conserved) because the global abundance 
of C is fixed, and the distribution of the grains has to com- 
pensate for the small number of VSGs with respect to the 
higher number of BGs. In Figs. ll7l and [T%l we show the effect 
of varying be on the emission of two young dusty SSPs of 
10 7 and 10 8 Myr both for the metallicities Z = 0.02 and 
Z = 0.008. At increasing be, the AIB emission increases, 
whereas the FIR emission goes in the opposite way. In our 
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Figure 17. PAH emission as a function of the parameter be- 
The MIR-FIR region is shown for two young dusty SSPs of 10 7 
and 10 8 Myr and metallicity Z = 0.02. Four values of be are 
taken into account, going from • 10 — 5 u ntil 3 ■ 10~ 5 , the upper 
limit f or the Ry = 5.5 extinction curve IWeingartn er &; Drainel 
l2001al) . The other para meters are kept fixed, R = 1, r = 35 and 
the ionization model of Wcingartncr & Drainc 1 2001b). 
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Figure 18. PAH emission in function of the parameter be- 
The MIR-FIR region is shown for two young dusty SSPs of 
10 7 and 10 s Myr and metallicity Z = 0.008. Three values of 
be are taken into account, going from ■ 10 — 5 until 2 ■ 10 — 5 , 
the uppe r limit for the average LMC e xtinction curve, as pro- 
posed in IWeingartner fc Drainel J2001al) . The other parameters 
are kept fixed, R = 1, r = 35 and ionization model of 
IWeingartner fc Drainel teOOlbl) . 



library, four values of be for solar or super-solar metallici- 
ties are considered, three values are taken into account for 
the LMC metallicity, and only one v alue is explored for the 
SMC a nd lower metallicities (see also lWeingartner fc Drainel 
12001 ah. 

PAH ionization: together with be, there is another pa- 
rameter affecting the emission of PAHs to be taken into ac- 
count, in order to match also the complicate profiles of the 
AIBs, i.e. the ionization state of the PAHs. The profile of the 
bands depends indeed on the PAH ionization, because the 
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Figure 19. PAH emission as a function of the ionization model. 
The MIR-FIR region is shown for two young dusty SSPs of 10 7 
and 10 s Myr and metallicity Z = 0.02. Three models of ionization 
are taken into account as indicated. 



optical properties are different for ionized and neutral PAHs 
as been discussed in Sect.|3] The ionization is calculated for 
each spherical shell of a MC. It strongly depends on the age 
of the SSP: passing from profiles of ionized PAHs for the 
youngest populations (because of the strong radiation field) 
to those of neutral PAHs at old ages (because of the weaker 
radiation field). 

In order to fully explore the space of parameters, 
three ionization model s are co nsidered. The first one is by 
IWeingartner fe Draind i2001fj) . The second model strictly 
follow s the same ionization profile calculated in lLi fc Draind 
fcOOlt) for the diffuse ISM of the MW. The third model sim- 
ply takes into account only the optical properties of neutral 
PAHs. In this way we can study a sequence of ionization 
states going from strongly ionized to almost neutral. For 
the second and the third model there is no change of the 
profile of the AIBs at varying the age of the underlying SSP 
as instead it happens in the first model. In Fig. 1191 we com- 
pare the results for the three ionization models showing the 
emission of two young dusty SSPs with different age, i.e. 10 7 
and 10 s Myr, and the same metallicity Z — 0.02 (all other 
parameters are kept fixed). The PAH ionization has no effect 
on the UV-optical and FIR emission because the distribu- 
tion of the grains is the same and their optical properties 
in these wavelength intervals do not depend on the ioniza- 
tion state. In the MIR we can notice that, as expected, the 
models with more ionization have more enhanced AIBs at 
6, 2 and 7.7 /xm and 8.6/xm, and weaker AIBs at 3.3 /im and 
11.5 /im. Clearly, the neutral model, has the strongest AIB 
at the 3.3 /im. 



8.2 Evaporation of the MCs 

Before using the SED of young dusty SSPs to derive the 
SED of galaxies, an important aspect has to be clarified. 
How long does it take to evaporate the MC surrounding a 
new generation of stars? The subject must be clarified in 
advance, because what we have been doing so far is only to 
take the energy emitted by a SSP and inject it into a MC 
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Figure 20. SEDs of a SSP of metallicity Z = 0.02 from the age of 
10 Myr (dashed line) to the age of 35 Myr (dot-dashed line). Two 
intermediate ages are shown of 20 Myr (continuous line) and 29.5 
Myr (dashed line), respectively. The parameter to is set equal to 
15 Myr. 
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Figure 21. The same as Fig.|3Jbut for the metallicity Z = 0.004. 



to be re-processed by dust: we have not specified the evo- 
lution of the MCs, which are most likely to diss olve by the 
energy input of the stars underneath. Following ISilva et alJ 
(1998), the process can be approximated by gradually low- 
ering as function of time the amount of SSP flux that is re- 
processed by dust and increasing the amount of energy that 
freely escapes from the MC. This is, of course, a rough ap- 
proximation of the real situation in which the energy input 
from supernovae and massive stars gradually sweep away 
the parental gas and dust. A more realistic approach should 
take into account the energy input coming from supernovae 
explosions and massive stars and compare it to the binding 
energy of the gas. 

Let us denote with /dssp the fraction of the SSP lu- 
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minosity that is reprocessed by dust and with to the time 
scale for a MC to evaporate 

C 1 t^t 
/dssp = < 2 - t/to t <t^ 2t (40) 

[ t > to 

Accordingly, the fraction of SSP luminosity that escapes 
without interacting with dust is fssp = 1 — /dssp- 
The parameter to will likely depend on the properties of the 
ISM and type of galaxy in turn. Plausibly, to will be of the 
order of the lifetime of massive stars, say from 100 to 10 
Mq, that first explode as supernovae and whose ages range 
from 3 to 100 Myr. However, the real timescale is likely to 
be much close to the shortest value, i.e. the lifetime of the 
most massive supernovae in the SSP, for a low density envi- 
ronment like spiral galaxies with a low star formation rate, 
while it will be much close to the longest value for a high 
density environment like the star forming central regions of 
starburst galaxies. For the sake of illustration, in Figs. 1201 
and 1211 we show the gradual transition of SSPs with differ- 
ent metallicity, Z — 0.02 and Z = 0.04, respectively, from 
fully embedded in their original MCs to free of dust. In all 
the cases, the optical depth is ry — 35 and R = 1, b c is fixed 
at the high value holding for MW a nd at the only available 
value fo r the SMC and, finally, the IWeingartner fc Draind 
ionization model is fully used. The SEDs are dis- 
played at four different ages, 10, 20, 29.5, and 35 Myr. The 
parameter to is fixed to to = 15 Myr. 

At the age of 10 Myr, the SSPs are fully enshrouded by 
dust and no evolution of the fluxes is still simulated (dotted 
lines in Figs.^Ul and|2TJ- At the age of 35 Myr, the SSPs are 
no longer surrounded by MCs, their SED is the one of a bare 
classical SSP without dust contribution of extinction and 
emission (dot-dashed lines in Figs. 1201 and l21H . For the other 
two values of the age the SED is intermediate to the previous 
ones. The trend is as follows: at the beginning the stars 
are heavily masked but when the age of the SSP enter the 
range between to ^ t ^ 2to part of the energy is transferred 
to the optical range and part remains in the IR. As the 
stars become older, the quantity of energy emitted in the IR 
diminishes and eventually disappears when t ^ 2to. 



9 COMPARISON WITH OBSERVATIONAL 
DATA 

To test our model for dusty SSP we compare the theoret- 
ical SEDs with the observational data for the central star 
forming regions of some local starburst galaxies. A sim- 
ple picture aimed at explaining most of the obse rvational 
constr aints for starbursters has been proposed bv ICalzettil 
feOOlf) . The newly born stars, likely still embedded in their 
parental MCs, are concentrated in the centre of the galaxy, 
surrounded by long-lived less massive stars wrapped by a 
diffuse dusty medium. Clearly, the real distribution of stars 
and dust will be much more complicated, with clumps and 
asymmetries. If the burst of star formation is strong enough, 
the flux coming from the central region of the galaxy (small- 
aperture data) is dominated by the obscured light coming 
from the new generations of stars. In this way, if we limit 
ourselves to study small central regions of starbursters, to 



a first approximation we can try to reproduce the observed 
SED with an obscured SSP of suitable mass. It is clear, 
however, that if the burst lasts lo nger than the life of the 
most massive stars (Calzetti 2001), the stellar populations 
in the central region under examination are no longer coeval. 
Rather they are a mix of spatially concentrated very young 
stars superimposed to older and more dispersed foreground 
stars. The longer or the weaker the burst of star formation, 
the higher is the influence of the older generations of stars on 
the resulting SED and the weaker the physical significance 
of the SSP approximation. Nevertheless, despite of these re- 
marks, the comparison of our SSP SEDs with the SED from 
the central dusty star forming regions of local starbursters 
is the best tool to our disposal to validate the model. 

9.1 Central region of Arp 220 

Arp 220 is a very br ight object in the Local Universe, clas- 
sified by ( 1966) as a peculiar galaxy. The interest in 
this galaxy became very strong when the IRAS satellite 
revealed that it is a powerful source of FIR emission. To- 
day, we know that Arp 220 is the nearest and most thor- 
oughly studied example of UltraLuminous InfraRed Galax- 
ies (ULI RG), objects with a n infrared luminosity Lir ^ 
10 12 L Q JSanders et al.lll988l) . The redshift of the galaxy, 
taken from the NED database is z = 0.01813, which for the 
Hubble constant Ho — 72 km/s/Mpc yields a distance of 
about 76 Mpc. Thi s value is fu l ly con sistent with the dis- 
tances proposed in ISoif er et al. l lll987fl an d in ISpoon et alJ 
J2004h . IWvnn- Williams fc Becklinl (Il993l) showed that al- 
most all the MIR flux of Arp220 comes from a small cen- 
tral region of about 5 aperture. This concentration of the 
MIR emission has been confirmed by ISoifer et al.l <ll999T) 
comparing the fluxes measured varying the beam size at 
a fixed MIR wavelength. We take the 5 central region of 
the galaxy, which c orresponds to a diameter of about 2 kpc 
jTakaei et alJl2003l) . to test our SSPs. 

In Fig. 1221 we show the results of our fit obtained us- 
ing dusty SSPs, with high optical depth tv = 35, compact 
structure with R — 0.5 and the SMC extinction curve. The 
low metallicit y extinction curve we have selected agrees with 
the results bv lTakagi et alJ i2003ft . Only with a mixture of 
grains, like the one of the SMC extinction curve, we were 
able to reproduce the deep absorption features of silicates, 
the low emission of PAHs and the optical/NIR continuum. 
So, the fit seems to require relative proportions of grains 
favouring the grains of silicates with respect to the graphitic 
ones, like in a flat SMC extinction curve, even if the average 
metallicity of the environment could be easily supersolar. 
To clarify the issue we should apply th e EPS to chemica l 
models in which the dust abundances toweklll998l [2005) 
are followed in detail. As this step is not yet accomplished, 
work is in progress. The best fit is found for the age of 30 
Myr, a mass of 0.4 • lO 11 A/0 and both solar and supersolar 
metallicity, thus suggesting the present of a high metallicity 
environment as shown in the right panel of Fig. 1221 

9.2 Central region of NGC 253 

NGC 253 is a bright starburst galaxy belonging to the Sculp- 
tor group, in the direction of the South Galactic Pole. It is 
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Figure 22. Left Panel: Best fit of the emission from the central region of Arp220. We show three SSP s with the ages of 10 , 30 and 
60 Myr (clotted, contin uous and dashed lines respectively). The best fit is for 30 M yr. Dat a are taken frompanders et alj ll988t B, g, r , 
n. lCarico et"aflll99(j J, H, K, L). The MIR fluxes a re fro mlSpinoglio et al.Hl995t MIRI. ISmith et ail I1989L MIR^. lKlaas et alj <1997t 
from MI R to FIR) . IWvnn- Williams fc Becklinl Jl993l M I Rt.lTran et alj fcOOll 5 - 16 Mm ISOCAM-C VF spectrum). FI R and radio data 
are fro m lRigopoulou et alj Il996l . sub-mml . lEales et all ll98St UKIRT s ub-mm), |Punne et alj J200ol. SCUBA sub-mm) , Ibunne fc Ealesl 
l200ll . SCUBA sub-mm). Other data are taken from Ispoon^taLl J2004I) , Right Panel: The same as in the left panel but for fixed age 
and varying metallicity. We show three metallicities: Z = 0.02, Z = 0.05 and Z = 0.1 (continuous, dashed and dotted lines, respectively). 



perhaps one of the nearest Galaxies besides those of the Lo- 
cal Group. The redshift is z = 0.000804 (NED), correspond- 
ing to the distance d — 3.35Mpc for Ho — 72 km/s/Mpc. 

The central region we have selected is within 15" aper- 
ture and its SED is centered around the peak of MIR 
emission. For this re gion there ar e avail able data on the 
PAH emission from ISturm et alj fcOOCt ISO-SWS) and 
IForster Schreiber et al.1 i2003L ISQCAM) that agree each 
other, thus confirming the accuracy of the calibrations and 
setting a strong constrain t on the global fit. Th ere are some 
problems with OAO data llCode fc Welcbll982tl covering the 
UV range, because of the large aperture 10 . We scaled the 
OAO data to the aperture size in usage here. It is clear that 
a precise correction would require the knowledge of the true 
UV profile of the galaxy, which is unknown. For this reason, 
even if the corrected UV level of flux is reproduced, we will 
give less weight to the UV region looking for the global fit. 

In the right panel of Fig. 1231 we show our best fit ob- 
tained with a SSP of solar metallicity, ry = 25, R — 1, 
MW extinction curve for dense regions (Rv = 5.5) and 
average-low value of the abundance of carbon in the very 
small carbonaceous grains (6c = 1.0 • 10~ 5 ). The ionization 
model which best agrees with the observations is the one 
with a dominant neutral component of PAHs. The mass of 
the population of stars is about 7 • 1O 8 M0. In the left panel 
we show an effect of a sort of age- mass degeneracy: if we can 
obtain a good fit for a given age and mass, an equally good 
fit is possible lowering both the age and the mass, because 
a younger population is more luminous. In spite of this un- 
certainty, there should be a young population of stars with 
age of about 30 Myr and total mass of about 5 ■ 1O 8 M0. 

9.3 Central region of M82 

The starburst galaxy M82 with its companion M81 was 
discovered in 1774 by Johann Bode and independently 



by Mechain in 1779. He reported them to Messier, who 
added these two objects to his catalog. The central re- 
gion of the galaxy, which l ikely suffered a strong inter - 
action about 10 8 years ago iForster Schreiber et al]l2003l) 
with its companion M81, shows a remarkable burst of 
star formation ac tivity. The classical distanc e to M82 is 
3.25 ± 0.20 Mpc jTammann fc Sandagdll968l) that makes 
of M82 the nearest and most stud ied starburst galaxy (see 
IShopbell fc Bland-Hawthornlll998l . and references therein). 

For M82 there are many available data for the PAH 
emission in the MIR region, to which we pay particular 
attention. They are very useful to constrain the fit to the 
observed SED of the central region. ISO-SWS data by 
ISturm et al.1 i200Ct private communication) cover a wide 
spectral region going from 2.4 to 45/im. Different SWS full 
grating scans that are observed with different aperture sizes, 
going fro m 14 x 20 to 20 x 33 are patched together. 
Recently, IForster Schreiber et al.l d2003l . private communi- 
cation) observed M82 in the MIR range between 5.0 and 
16fim with ISOCAM on board of ISO. The total field of 
view was 96 x 96 . We have chosen to examine the central 
M82 region of 30 x 30 . For this coverage we can use the 
IForster Schreib er et alj j20 p3l) dat a (Core plus Disk obser- 
vations), while Istorrnetal T T200CI) data with smaller aper- 
tures have been properly corrected. 

In Fig. 1241 we show our best fit, obtained with ry = 
10, R = 1, metallicity Z = 0.008 and the LMC extinction 
curve. The carbon abundance in small grains is equal to the 
average be = 1.0 ■ 10 _J of current estimates for LMC. The 
best ionization state of PAHs is the average ionization model 
for the MW diffuse ISM. All the spectral regions are well 
reproduced, the major uncertainties being the 18fim silicates 
absorption feature and the far UV region where we have no 
data. The mass obtained from the fit is M = 8 • 1O 8 M0. 
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Figure 23. Left Panel: Best fit of the emission from the central region of NGC253. We show three SSPs with ages of 10, 40 and 95 Myr 
(dotted, continuous a nd da shed lines, respe c tively ). The best fi t is for 40 My r . Data for the central regi on are taken and adapt ed from 
Code fc Welch! <1982t UV). ISpinoglio et alj Jl995l J, H, K, L). Ijarrett et alj J2000j. J. H. K - 2MASS ), iMoshir fc et"aD Il990t IRAS), 
Chini et alj <1984l sub-mmi . ISturm et alj feOOCt private communication), iForster Schreiber et ail J2003I. private communication). Right 
Panel: The same as in the left panel but at varying both the age and the mass of the galaxy for the same metallicity. 
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Figure 24. Modelling the SED of the 30" X 30" central re- 
gion of M82. We show three SSPs with ages of 10, 25 and 
60 Myr (dashed, continuous and dotted lines, respectively). 
The best fit is for 25 Myr. Data for the central region of 
this g alaxy are taken from: NED online database, |jcjmson| 
11963. UBVRD.Ilchikawa et all Il995l JHK) . ISoifer et alj 119871 
IRASI.lGolombek et alj Il988l. IRAS ). iRice et alj <198sT lRAS). 
iKlein et alj Il988l FIRand sub-mm).ISturm et alj I2OOOI private 
communication), Forstcr Schrei ber et al .1 i2003l private commu- 
nication) . 
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Figure 25. Modelling the SED of the 15" X 15" central region of 
NGC1808. We show three SSPs with ages of 30, 95 and 250 Myr 
(dashed, thick continuous and dotted lines, respectively). The the 
best fit is for 95 Myr for be = 1.0 TO 5 . Another good fit (thin con- 
tinuous line) with different be is also shown. D ata for this g alaxy 
are ta ken in the NED onl i ne dat abase fr om i Gricrsmit h et alJ 
Jl98l JHK) Ijarrett et alj feOOOt JHK ). iGlasj Jl973t HKL) 
i Frogel et alTjl982l. MIR. 8-13 urn) (MIRlForster Schreiber et all 
l2003T ~p"rivate communication) and lMoshir^^t^n * i^990r TRAS). 



9.4 Central region of NGC 1808 

NGC 1808 is a starburst spiral galaxy in the southern hemi- 
sphere sky. The redshift of the galaxy is z = 0.003319 (NED) 
which for Ho = 72 km/s/Mpc corresponds to a distance of 
about 13.82 Mpc. 

For NGC 1808 we have chosen the region 15 x 15 in 
the centre of the galaxy. In Fig. 1251 we show four models, 
all calculated for the metallicity Z — 0.008, optical depth 
tv ~ S, scale radius R — 1. The best fit is obtained for the 



age of 95 Myr and a carbon abundance in the log-normal 
populations of VSGs of be = TO ■ 10 -5 . A good fit is also 
possible for a lower abundance be — and slightly different 
age of 85 Myr (see Fig. 1251 . In all models, the total stellar 
mass is M = 5 • 10 9 M Q . The case of NGC 1808 is more 
problematic than the other galaxies we have examined in- 
sofar. First, there are no data for the optical region with 
the desired aperture and the aperture-correction applied to 
other optical data published in literature is too problematic. 
Second it is difficult to simultaneously reproduce both the 
flux level of the PAH emission and the absorption features 
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of silicates at 9.7/xm and 18/im. In our fits the absorption 
features are too deep. The ideal solution would require an 
obscuration curve with different abundance of silicates. Fi- 
nally, some data at sub-mm wavelength would be needed to 
better constrain the fit. Anyway, the results we have derived 
suggest a population with an average age of about 100 Myr. 
In the case of this galaxy the SSP approximation is perhaps 
too crude and it should be more appropriate to use of a mix 
of stellar populations, thus perhaps getting a better fit of 
both the absorption and emission features in the MIR. 



10 DISCUSSION AND CONCLUSIONS 

In this study we have presented detailed models of extinction 
and emission of the diffuse ISM and the application of the 
most suitable one to the calculation of a library of young 
dusty SSPs embedded in their parental MCs. 

We have analyzed the properties of a dusty ISM hav- 
ing as guide-line of the strategy the cross-checking of the re- 
sults for extinction and emission and iterating the procedure 
to contrive the many, somehow unavoidable parameters at 
work, so that not only unrealistic solutions have been ruled 
out (they could by generated by the sole fit of the extinction 
curve) but also additional information on the whole problem 
has been acquired. 

We have adopted two possible models for the dis- 
tribution of the grains, shortly named MRN and WEI, 
and explore the consequences. In brief, the MRN-model 
is a multi-power-law, whose number of components varies 
with the grain type, whereas the WEI-model is a com- 
plicated analytical law suitably corrected for the smaller 
carbonaceous grains by the sum of two log-normal terms. 
Again, we followed two methods to derive the tempera- 
ture distribution of VSGs. The first method is a combi- 
nation and modification of wh at has been proposed from 
iGuhathakurta fe ^Jraine] (^)89|) fo r silicate and graphite 
grains and from 1pTge^etalTiu98^ for PAHs, in which we 
introduc ed the calculation of the ioniza tion state of PAHs 
following lWeingartner fc Drainei ll2001b|) . It has been shortly 
indicated as the G DP-mode l. The second method, due to 
iDraine fc il (1200 ll) and lLi fc Drainej fcOOlll . is the state-of- 
the-art of these problems as it stands on the exact statistical 
solution of the problem. We have referred to it as the LID 
model. Using the LID and/or GDP emission model coupled 
with MRN and/or WEI model for the size distribution we 
have calculated the extinction and emission properties of 
the average diffuse ISM of the MW, LMC and SMC and 
compared it with observational data. 

Looking at the results from GDP+MRN and 
GDP+WEI and their comparison with the LID+WEI, we 
adopt the GDP+WEI model as the most suitable descrip- 
tion of the properties of the dusty ISM to be used in dusty 
SSPs and in the model of galax y SEDs in presence of inter- 
stellar dust JPiovan et al.ll2005fl . 

We have applied the model of dusty ISM to calculate 
the SEDs of young SSPs, in which the effects of dust are 
taken into account. These are particularly relevant in mas- 
sive, young stars when still embedded in their parental MCs. 
We propose a description in which stars, gas and dust in a 
MC obey three different density profiles each of these charac- 
terized by its own scale length. In this context, we derive the 



optical depth of the MC and apply the ray tracing method 
to solve the problem of radiative transfer across the MC. We 
have also explored the parameter space of the MC models. 
Four parameters have been taken into account, namely the 
scale radius of the MC and its optical depth (the geometri- 
cal parameters) , the abundance of carbon in the log-normal 
populations of very small carbonaceous grains and the ion- 
ization state of PAHs, both strongly influencing the PAH 
emission. Basing on this, we have generated an extended li- 
brary of young dusty SSPs of different chemical composition 
and age specifically taylored for galactic population synthe- 
sis, the final goal of our study. We tested our SEDs trying to 
reproduce the observational data of star forming regions in 
the centre of local starburst galaxies, reaching a good agree- 
ment between theory and observation. Libraries of SSPs at 
varying these parameters are made available to public for 
future use 2 . 
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